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ABSTRACT 


The trade-off between false alarm and detection probability is a 
fundamental challenge in the automated detection of contacts in passive sonar 
systems. A common approach is the application of high-gain processing 
followed by successive classification criteria. Most classification schemes (e.g., 
matching of signature) are complex and tailored to specific target types. By 
contrast, the Doppler effect is readily observed in all contacts with discrete tonals 
and relative velocity to the receiver. This thesis demonstrates that the Doppler 
effect can be exploited to improve the detection process by filtering out contacts 
that do not exhibit these characteristics. Cross-correlation (matched filtering) of 
contact LOFARgrams with templates generated jin situ with Doppler 
compression/dilation is used to achieve this. Velocity information and the range 
of the Closest Point of Approach (CPA) are estimated in the process. A 
detection algorithm is developed in MATLAB, and the radiated acoustic 
signatures of overflying airplanes are recorded at a ground station. In the 
analysis of six propeller and four jet airplanes, the program successfully identifies 
the passage of all six propeller airplanes with four incidences of false alarm, due 
in one case to a jet airplane. Velocity and range estimates are also within 


expected values. 
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I. INTRODUCTION 


A. BACKGROUND 


The Seaweb program uses a network of underwater acoustic modems to 
communicate without the use of cables. These networks, once laid, provide 
round-the-clock presence in the area of interest. Recent work explores the 
possibility of using these underwater modems for a surveillance role in shallow 
waters and harbor entries. To achieve this, the components used in the existing 
network may need to be supplemented with additional hardware and software to 
sense the passive acoustic signature generated by maritime vessels. In addition 
to the provision of continuous monitoring of the environment around the modem 
for passing vessels, the collected acoustic signatures or shorter contact reports 
can also be transmitted ashore for analysis. 

In a multiple-sensor network, the amount of information that needs to be 
processed by the onshore operator increases with the number of sensors. Even 
before analysis of any collected signature can be done, the validity of a contact 
as a maritime vessel, as opposed to other competing noise sources, has to be 
determined. Automation of detection is the obvious solution in a contact-rich 
environment. However, the fundamental trade-off between the probability of 
detection and the probability of false alarm remains the challenge in the design of 
the automation logic and system. Setting too low a threshold results in a high 
false-alarm rate, while the converse results in a high chance of missed detection 
of an otherwise valid contact. 

The traditional detection system trigger is commonly based on the 
received signal being above a predetermined threshold, which is usually 
determined by operational requirements. The selection of an appropriate 
threshold that does not result in too many false alarms then becomes the 
challenge. Techniques such as the Neyman-Pearson criterion may be used to 
compute the threshold for a desired false-alarm rate (Ziomek, 1994). 
Alternatively, the use of a second independent property of the contact of interest 
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can serve as a second-level filter to weed out false alarms before declaring a 
detection decision. This thesis explores the use of the Doppler effect, which is 
visible in moving contacts with narrowband frequency components, as a second- 
level filter to decrease the probability of false alarm in the automated triggering 
mechanism. In addition, the velocity and passing range, which are important 
contact parameters, are estimated from the Doppler shift exhibited by the moving 
contact. 


B. RESEARCH MOTIVATION 


The trade-off between the false alarm and detection probability is 
fundamental in radar and sonar. (Chevalier, 2002) 


A common approach to the detection of passive contacts is to apply high- 
gain processing to the raw signal to increase the signal-to-noise (SNR) ratio. 
The processed signal is then compared to a threshold value. 

Statistically, if the difference between the contact’s mean signal level and 
that of the background noise is big, the probability of false alarm (Pi) can be 
reduced by setting the detection threshold higher without significantly increasing 
the probability of non-detection (Pra) of a true contact. However, in cases where 
this difference is small (as in case of a weak contact), the level at which the 
threshold is set involves in a trade-off between reducing Py, and Png. 

Many of the currently used classification methods are complex and mostly 
aimed at specific contacts (e.g., matching of contact signature). While the use of 
the Doppler shift is well established in the fields of active radar and active sonar 
for deriving velocity information on the contact, it is not as straightforward in the 
passive case due to the lack of a known transmitted waveform with which to 
cross-correlate the return echo. However, the Doppler effect is still an attractive 
classification criterion for the detection of moving contacts as it is prevalent in all 
moving contacts showing tonal lines. The objective of this research is to 
demonstrate the possibility of using the Doppler effect to discriminate legitimate 


moving contacts from background clutter, which will in turn allow the threshold for 
auto detection to be set closer to the noise floor without significantly increasing 


the false-alarm rate. 


C. RESEARCH SUMMARY 


This thesis is limited to the demonstration that the Doppler effect from 
moving contacts can be used to discriminate a moving contact with distinct tonals 
from competing sources of noise. The velocity and range of the contact are also 
estimated to a limited degree of accuracy. 

The theoretical framework of the tactical trade-off between probability of 
detection and probability of false alarm is first developed. The physics of the 
Doppler effect are then presented. Incorporation of the Doppler effect as a 
second-level filter is demonstrated as a means for lowering the detection 
threshold without increasing the probability of false alarm. 

The Doppler effect is commonly and routinely used in the time-frequency 
analysis of contacts in passive sonar. It is apparent on a spectrogram and is 
readily identifiable by a trained operator. As the goal of this research is to 
automate the process, the challenge is to develop a signal processing algorithm 
that not only recognizes the Doppler effect when it is present, but also analyzes 
the Doppler curve as a trained operator would do. Once this algorithm is 
programmed, a suitable candidate target has to be selected to test the fidelity of 
the algorithm. 


Acoustic data of overflying airplanes on approach for landing at the Monterey 
Peninsula Airport were collected at a land station against a background 
consisting of natural and manmade noise, such as automobile traffic. An 
algorithm was developed and programmed in MATLAB to search the collected 
data and detect a passing airplane by identifying the Doppler effect. The results 
were compared against airport records to check the fidelity of the program. The 
algorithm is shown to be extremely successful in discriminating between loud 
contacts which do not exhibit an easily discerned Doppler shift (jet airplanes) vis- 
a-vis those exhibiting a clear Doppler shift (propeller planes). 

Finally, recommendations for further work to apply the algorithm to an 


underwater environment and passive sensing systems are presented. 


ll. THEORY 


A. THE TACTICAL TRADE-OFF BETWEEN PROBABILITY OF 
DETECTION AND PROBABILITY OF FALSE ALARM - A STATISTICAL 
ANALYSIS 


1. Detection Performance 


The detection process involves a decision about the presence of a contact 
through the processing of its signal (acoustic, radar returns, etc.) against the 
ambient noise of the background medium. A detection decision is declared if the 
received signal level exceeds a predetermined threshold. In the sonar context 
where logarithmic quantities are used, the detection decision can be represented 


by the following condition in the sonar equation: 


Signal — Noise + Gain > Threshold. Eq. 2.1 


This “signal” is not limited to a time domain or frequency domain signal but 


can be extended to mean any measurable output of the sensor system. 
The performance of a system is characterized by two statistics, namely: 


Probability of Detection, Pg, defined as the probability that a processed 
signal exceeds the predetermined threshold when the contact of interest is 


present. 


Probability of False Alarm, Pr,, defined as the probability that a processed 
signal exceeds the predetermined threshold even though no contact of interest is 


actually present. 


It is immediately obvious that the correct selection of the predetermined 
threshold is the key to the system’s performance. A conservative decision 
criterion designed to minimize false alarms translates to setting a high threshold. 
This means that only strong signals result in a positive detection call, giving a 
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small Py, and a correspondingly lower Pg. On the other hand, an aggressive 
criterion designed to maximize the probability of detection would use a low 
threshold to increase Py at the expense of a high Pr. In Figure 1, the “signal” 
being evaluated is the power output as a function of time. A lower threshold 
increases the number of false alarms triggered by noise spikes, while a high 
threshold reduces the number of false alarms by noise spikes but also results in 


legitimate signals being missed. 





False alarm 
oe Detected 
QY 
_> | Undetected signal 
) . 
= signal ms 
a) 
By POO eS ena po = — hreshold 
so) 
iD) 
fs 
o 
9 Mean noise 
as level 
Noise 
Time 
Figure 1. __ Illustration of Effects of Varying Threshold on Detection and False 


Alarm (From: Wagner, Mylander, & Sanders, 1999) 


2. Optimizing Performance 


In both sonar and radar applications, optimization of system performance 
is done such that Py is maximized for a given upper limit of Pf. The aim is to find 
a balanced and acceptable threshold level (false-alarm rate). This decision is 
influenced by, among other factors, the operational imperatives of the system 
favoring one criterion or another. In some cases (e.g., detection of enemy 
submarines), a non-detection might have fatal consequences, whereas false 
alarms have less serious consequences, such as a waste of time and effort. 
Here, an aggressive decision criterion designed to maximize the probability of 


detection would be favored. 


To examine the trade-off between the probability of detection and the 
probability of false alarm in more detail, consider a case where the received 
signal is Gaussian regardless of whether a contact is present or not. The actual 
statistics of the received signal in both cases depend on the nature of the signal 
of interest, the noise, and the way in which the raw signal is measured. 
However, any assumed statistics serve to illustrate the trade-off. In the 
discussion which follows, let s,(t) be the amplitude of the received signal. The 
signal due to the contact of interest is s(t), and the selected detection threshold 
against which the signal is compared for a detection decision is 7. Assume 
further that the noise-only signal is stationary and given by n(t) which has a 


mean of uw, and a variance o*. If the contact is present, assume that the 
received signal has a mean of w,,, and a variance of o°. The following two 


Cases May OCCUr. 


a. Contact is Absent 


In the case where there is no contact of interest in the environment, 
the only signal received by the sensor is Gaussian noise, n(t). Hence, s,x(t) = n(t) 
where s;(t) ~ N ( “iv, o” ). Since there is no contact, any detection decision is a 


false alarm. The following equation applies and is illustrated in Figure 2. 


Py =Pr 8, )27 |8,(0)~N (Uyio*) } 


1 . 2 (X — Hw y 
oN2n peep 20° ve 


1 (7 — He) Eq 2.2 
oO 





where © is the cumulative normal function. 
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Figure 2. Illustration of the Noise-Only Case 


b. Contact is Present 


When a contact of interest is present, the received signal is the sum 
of the signal from the contact s(t) and the noise n(t). The probability of detection 
is thus the probability of the sum of signal and noise exceeding the threshold, 


where 


Sn(t) = s(t) + n(t) and s,(t) ~N (nes, 07 ) 
ae =Pr Ss, ¢)27 |s,¢€)~N iceO) } 
Hnss j 
= app fexe CMSs oe 


-1-0( To #us | Eq 2.3 


Oo 





This is illustrated in the following figure. 


Target Present yy 
Received Signal + U 
Noise 

Sx (t) ~ N ( LNsS , o*) / > 
Thresnotl 


Figure 3. __ Illustration of Target+Noise Case 


Putting both the above cases together, the following illustration, 
with both of the above probability distributions, shows the situation encountered 


by the sensor in the presence of noise. There will inevitably be false alarms 


caused by noise only regardless of the selected threshold and the link between 
Pg and Pr, can clearly be seen. 


Target absent 
Noise only 


Sx (t) ~ N (Law 0) 


Target present 
Signal + Noise 
Six (t) ~ N (Lines ,O) 







Py 





I 
Ln LIN+S 


Figure 4. | Sensor Detection in Background with Noise 


If the “costs” associated with false alarms and missed detection can 
be quantified and are known, an optimized threshold can be computed 
mathematically (Washburn, 1996). 


Summarizing all possible cases results in the flowing matrix: 











Target Present Target Absent 
Sensor indicates Pra (False Alarm / 
Py (True Positive) 7” 
Target Present False Positive) 
— 1- Py 
Sensor indicates 1- Pra 
(Missed Detection / 
Target Absent (True Negative) 
False Negative) 

















Table 1. Target Matrix 
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3. Receiver Operating Characteristics (ROC) Curves 


The ROC curve is a plot of Py against P;, for a sensor system. It was first 
used during World War II for analysis of radar signals before being employed in 
signal detection theory (Green and Swets, 1966). These curves show the trade- 
off between Py, and Py, for different signals of interest and noise statistics as the 


detection threshold is changed. 





From Equations 2.2 and 2.3, P, -1-0[ Sf) and P, -1-0[ 24s) 
os o 


Let the normalized threshold be vole The equations for the probability of 
oO 


detection and false alarm can then be written p,=1-0[ x | and 
oO 


rae -1-0[x-41). This makes it clear that the difference in the probability of 
oO 


detection and the probability of false alarm depends on the difference in the 


values of Een and Py as well as on the choice of threshold. The detection 
oO oO 
_ 2 
index is defined in terms of the difference of these values as d -( ferme) . It 
Oo 


is a measure of the signal-to-noise ratio. Higher values of the detection index 
yield better detection performance for a given threshold value. At the same time, 
for a given detection index, Py and Py, increase with lower threshold values. 


These relationships are shown in the ROC curve below. 
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Figure 5. | Generic ROC Curves 


B. THE DOPPLER EFFECT 

The Doppler effect refers to a shift of the apparent signal frequency 
caused by the relative motion of the source and receiver. Valuable information 
about the contact’s trajectory can be derived from its Doppler shift. 
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1. Near the Closest Point of Approach (CPA) 





" oh 
! eX Velocity = v 
Ro: Range atCPA | orp 
mo) 
Receiver oO 


Figure 6. Passage of Airplane through CPA 


Consider the above figure, illustrating an airplane flying by a stationary 
receiver with constant velocity, v. If the time of CPA is taken to be zero, the 
component of the plane’s velocity towards the receiver along the line of sight at 


time tis given by: 


vt vt 


Vos = VSin 8 = yp—- = ———————— 
LOS r(t) fon? +R, EQ 2.4 


The apparent frequency received at the receiver at any time, 





Vios Vios 
pe co) —-EQ25 


qi Vios a= Vios vies Vios ) in (Yzos y? 
Cc Cc Cc Cc 


where fo is the frequency received at CPA when the plane has zero velocity 
along the line of sight and c is the speed of sound in the medium. 
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For c >> Vios , equation 2.5 simplifies to the following : 


f SF 1 4 eae EQ 2.6 


Rearranging and differentiating with respect to t yields : 


@Vios _ a (as EQ 2.7 
dt T% 


When the plane is close to CPA, O << 1, hence 


Vv 
Vios =VSin = vO = v tan 8 = —v— EQ 2.8 
R 


O 


a 


where ft is negative before CPA and positive after. 


Substitution of Eq 2.8 into Eq 2.7 yields : 


BWih 2... VO se oo . 
dt  — R are 


df _ WO he EQ 2.9 
dt R C 





The rate of change of frequency as the contact passes through CPA is 
proportional to its speed and inversely proportional to distance. Therefore, a 
measurement of the rate of change of frequency near CPA can be used to 


estimate the range at CPA. 
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2: Away From the CPA 


The received frequency of the contact approaches an asymptotic value 
when it is far away; hence, the center frequency can be computed from the 
opening and closing asymptotic values: 


f = by fi EQ 2.10 


+ 
oO 2 ) 
where fu is the closing (upper) frequency asymptotic value and fL is the opening 
(lower) frequency. The Doppler shift is then given by: 


ie ae on ea oe EG i 
Cc 


Rearranging, the speed of the contact, v, gives: 


cAf EQ 2.12 
Fs 


Through the measurement of fy and f,, the velocity can be determined, 
which can then be substituted back into EQ 2.9 to determine the CPA range Ro. 


3. Doppler — Range Ambiguity 


The recorded pressure vs. time signal in the time domain is a noisy 
waveform that conveys little information to the human eye and yields little 
information for easy analysis. Therefore, a Discrete Fourier Transform (DFT) is 
used to transform it into its corresponding representation in the frequency domain 
so that the amplitude of the component sinusoids can be extracted. However, 
the two parameters, time and frequency, are conjugated by the process of 


Fourier transform. 


ie: 


The process of transformation involves first multiplying the signal by a 
smooth curve (typically a Hamming window), following which, the DFT is taken. 
However, if the DFT for the entire sample is taken as one long data sequence, 
the variance of the noise can be large, making it difficult to observe the frequency 
components of the signal of interest. It is common practice instead to break the 
input signal into many smaller segments. Each of these segments is multiplied 
by the Hamming window, run through the DFT, and the resulting frequency 
spectra are then averaged to form a single frequency spectrum. This process of 
averaging reduces the noise variance, allowing interesting features of the signal 
to become prominent. However, the length of the DFT translates to frequency 
resolution. By taking segmented parts of the original signal, we trade frequency 
resolution to provide more “averages.” For a finite signal, the number of 
averages is governed in turn by the length (integration time) of each segment. 
Hence, if we desire a spectrum of high-frequency resolution, longer segments 
(longer integration time) would be required and vice versa. This results in a 


Frequency (Doppler) - Time (Range) ambiguity that is inherent in our analysis. 


C. CROSS-CORRELATION AND IMAGE DETECTION 


The concept of correlation involves the mathematical operation of two 


signals (received signal and template waveform) to generate a third signal. 
f(itweg(t)= [  F*(f wttte et coz 


When the two signals are different, the process is known as cross- 
correlation. The following is a simple illustration of the cross-correlation process. 


For two finite length real discrete sequences, x(n) and y(n), the cross-correlation 


is defined as : 
: N -1 
Rgy (in aes x(n y(n +m ») EQ 2.14 
n=0 
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x(n) is a received signal within which we wish to look for a certain 
waveform, y(n), the template waveform. This is achieved by taking the sum of 
the product of x(n) and y(n+m) for different values of m to give a corresponding 
value of Ry. By "shifting” y(n) by delay m, the amplitude of R,y at each "sampling 
position” is therefore a relative measure of how closely a particular segment of 
the received waveform x(n) resembles the template waveform y(n). When 
plotted, a peak occurs in Ryy if the template waveform is present in the received 
signal x(n). Extending the same principle, correlation can be used to look for a 
known pattern in an image. The pattern is correlated to different segments of the 
main image or different sample images are correlated to the main image to find 


the location of the best match or a sample of the best match. 


Cross-correlation involving a template signal is also known as matched 
filtering. 
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x(n) 





Figure 7. Illustration on Concept of Cross-correlation (After: Smith, 1997) 


lil. EXPERIMENTAL AND COMPUTER PROGRAM 


A. DATA COLLECTION 
1. Choice of Data Source 


For consistency with the assumptions made in the development of the 
equations introduced in Chapter Il, the contact should have a constant course 
and speed as it passes CPA. To establish ground truth, the contact should also 
be identifiable. In addition, the contact should also have tonals with Doppler 
shifts in the mid-frequency range, where they are best observable. For an 
assessment of the accuracy of the program’s estimate of velocity and range, the 
velocity and CPA should also be known. 


It is fortuitous that the Naval Postgraduate School (NPS) is located near 
the flight path of the Monterey Peninsula Airport. This makes it convenient to 
make recordings of airplanes passing overhead on their final approach to land. 
These airplanes meet the above criteria as they make minimal course and speed 
changes during this phase of flight. Distinct tonals from the propeller airplanes 
are present in the 100-1000 Hertz (Hz) frequency range. As will be shown in the 
next paragraph, the range at CPA can be estimated. The airplane’s velocity can 
be deduced based on the identity and type of airplane. These conditions made 
them good candidate targets for this study. 


2. Location of Recording 


An hour of recording was conducted on 07 Aug 2008 from the top of the 
tower of Herrmann Hall, NPS. The airplanes were either flying the final approach 
of the Instrument Landing System (ILS) glide slope of the Monterey Peninsula 
Airport's runway 10R (depicted by dotted lines in the aerial photo) or on the 
visual approach to the parallel runway 10L. 
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Figure 8. ‘Aerial Photo of the Naval Postgraduate School and Surroundings (From: Mapquest, n.d.) 
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R/W 10L 


R/W 10R 





Figure 9. Schematic Diagram of Target Track Relative to Fixed Observing 
Station 


Based on measurements taken from the above aerial photo, the distance 
between the tower of Herrmann Hall and the runway is estimated to be 1.5 
kilometers (km). For airplanes flying overhead on the 3° ILS glide slope, a 
passing distance of 80 meters (m) between the receiver microphone and airplane 
at CPA is approximated. Airplanes on approach to the parallel runway 10L are 
estimated to have a passing distance of 200 m. 
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3. Equipment and Setup 


Microphone power supply 


Microphone & Pre-amp - 
#8 


2 





Figure 10. Photograph of Equipment Setup used in Recording 





Microphone and Pre-Amp 


NI 9233 NI USB- 
DAQ 9126 Signal 
Converter 


Figure 11. Graphical Representation of Equipment Setup 


22 


The equipment used for the recording is listed as follows. The technical 
specifications are included in Appendix A. 


a. ACO Pacific Microphone. — Cartridge model 7046 & 
4012 —- : inch Preamplifier 


b. National Instrument NI9233 Dynamic Signal 
Acquisition modules with USB carrier. 


C. Dell D820 Laptop 


4. Processing of Collected Data 


Acoustic data are processed and stored in the National Instruments’ (NI) 
proprietary Technical Data Management - Streaming (.tdms) file format using the 
included NI LABVIEW SIGNAL EXPRESS version 2.5 software. This is 
subsequently converted and stored as an ASCII text file (.txt) using LABVIEW 
Express’ built-in file converter prior to the commencement of the post-analysis 


phase. 


B. MATLAB COMPUTER PROGRAM —- AUTO DOPPLER DETECTOR.M 
1. Program Overview 


MATLAB Program ‘auto_Doppler_detector.m’ was written to analyse the 
recorded data. The purpose of the program is to detect airplane passages and to 
determine if they are propeller airplanes through analysis of the detected energy 
level in conjunction with the presence of Doppler shift. The collected raw signal 
(uncalibrated sound pressure) is input as a .txt file. In addition to flagging the 
propeller airplanes that exhibit Doppler shifts, the program also returns the time 
of CPA, frequency of the strongest tonal, and an estimate of both the velocity and 
range at CPA. On completion of analysis, plots containing pertinent information 
on the passage and a summary of the detected peaks are produced to form the 


post-analysis report. 
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2. Program Structure 


The structure of the program is represented in the following flow chart. 


Measured 
signal input as 
txt file 


Calculate signal 
power spectrum 
density 


Summation of 
power in time 
period 













Is power > 
threshold? 





Yes 


Locate time of 
CPA and 


strongest tonal 





Calculate df/dt 


Generate template 
based on df/dt, fo and 
possible vel. 











Template 
matched Signal? 


Record and 
return velocity 


of best match 





Figure 12. Flow Chart of Program auto_Doppler_detector.m 


The power spectral density of the time domain input signal is estimated 


using the MATLAB built-in function ‘pwelch.m’ via Welch's method. By summing 


the total energy in the time bands, the program detects the presence of an 
airplane passage or increase in sound energy and flags the corresponding time 


period for analysis. 


A search is conducted within the flagged time period for the time with the 
highest energy. This time is associated with the time of CPA, tcp,, based on the 
assumption that the energy of a passing airplane is highest at CPA. 


The strongest tonal at time tcp, is identified and its frequency, fo is noted. 
This tonal forms the basis for Doppler analysis. The corresponding strongest 
tonals in the preceding time (before tcp,) and proceeding time (after tcp,) within a 
frequency span of Af are likewise determined. Af is calculated using Eq 2.11, 
using the highest expected contact speed. Through the tracking of these tonals 
during the passage of the target near CPA, the rate of change of frequency with 


respect to time is computed. 


Cross-correlation of the recorded signal’s tonal in the region of CPA 
against synthetically generated Doppler profiles of a range of hypothetical 
velocities with the determined df/dt is then carried out. The result of the cross 
correlation factor, Ry, at various velocities is tabulated. A high Ry, value is an 
indication of a good match and can thus be used to deduce the presence of the 
Doppler curve. CPA range is then computed based on the velocity of best 
match. 


3. Considerations in Design of the Program 


As the program is envisaged to reside in an underwater modem, the size 
of the program and the speed in which the results are obtained are 
considerations in its design. This led to the rejection of using pre-generated 
template signals of all possible contact velocities and passing ranges stored for 
correlation against the recorded signal. The library of templates covering all 
possible scenarios would require large storage space and accessing them would 
slow down the program due to their sheer number and time required to load 
them. Instead, to achieve a small memory requirement, templates are generated 
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in situ based on the measured frequency rate change and fo. This removes the 
need for any stored templates and greatly reduces data processing time. The 
MATLAB function ‘generate_template.m’ generates the required templates. This 
program, as well as all others written to support this thesis are provided in the 


appendices. 


4. Underlying Assumptions and Significant Decisions 


The accuracy, fidelity and robustness of any algorithm hinges on the 
assumptions made and the parameter estimate at key decision points. The 


underlying assumptions and their implications are discussed as follows: 


a. Estimate of df/dt Near CPA 


In order to keep the assumption in Eq 2.8 valid, (O << 1), the 
number of differences that can be taken to provide an estimate of df/dt in the 
vicinity of CPA was limited to one sample point before and after CPA in this 
study, i.e., di/dt is an average of only two estimates. This is due to the close 
range, Ro, of 80 m between the airplane and receiver microphone, high speed of 
the planes - between 30 to 70 metres per second (m/s) — and uncertainty in the 
estimate of CPA time tg (1 second) as a result of conducting the analysis in the 


frequency domain. This calculation is illustrated in Figure 13. 


44m 


For sin © = © = tan © (accurate to 1* decimal), © max * 
28.65°. For an airplane travelling at 50m/s, to keep O < 
80m © max, the period of ‘valid time’ is 0.88 second before 


and after CPA. 


Figure 13. Calculation of “Maximum Valid Time Period” 
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b. Tracking the Doppler Shift and Calculation of df/dt 


To calculate the rate of change of frequency in the vicinity of CPA, 
the program “looks” for the strongest tonal within a frequency range in the vicinity 
of fo in the time periods adjacent to CPA. The assumption that these strongest 
tonals are in fact due to the Doppler shift of fo subjects the analysis to possible 
error. Spurious noise or other signals from the contact with a higher amplitude 
than fo could lead to the misidentification of the tonal line and an error in the 
estimation of df/dt. This “risk” is mitigated in the program by limiting the 
frequency range over which the search is conducted. 


Cc. Flagging Inadmissible df/dt 


The sign of df/dt for a true Doppler shift will always be negative as 
the contact proceeds through CPA. This property is used as a validity check on 
the contact. In instances where a non-negative df/dt is encountered, the tonal is 


flagged and no further analysis is performed. 


d. Tonal Frequency Range 


From the analysis of the power spectrum of background noise and 
typical airplane passages, it is apparent that the level of background noise is 
considerably higher at the lower frequencies and falls off rapidly as we progress 
to the higher frequencies. It is also observed that Doppler shifts are best seen 
for tonals within the 100 Hz to 2000 Hz frequency range. In the analysis 
conducted, the search for tonals is limited to above 90 Hz. This is done to 
exclude possible interfering sources in the low frequency range where the 
background noise contains a strong, random broadband component which 
obscures the target. This lower frequency limit can be adjusted by changing the 


associated parameter in the program. 
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IV. DATA ANALYSIS 


A. OVERVIEW 


Data from a total of ten airplanes (both jet and propellers) were recorded 
during a one-hour period. The exact time of CPA was recorded manually as the 
airplane passed overhead. Identification of the target type was based on 
listening to the air-traffic control radio communications and subsequent 
verification with airport records (Flightaware, n.d). Post-processing of data 
included the audio playback of the events, analysis of pressure vs. time (time- 
domain analysis), and study of the targets’ tonals on the LOFARgram (Low 
Frequency Analysis and Recording gram, i.e., time vs frequency). Table 2 
summarizes the target events. This ground-truth forms the basis for evaluating 
the output results obtained from running the MATLAB program at different 
threshold settings. Time references to the passage of the airplanes are made in 
both real time and relative time measured from the start of the recording in 


seconds. 
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Relative Time 


















































Time of Event (s) Type of Target Airplane ID Remarks 
Ss 
Single 
1021 74 Cessna 182P N6800M 
propeller 
Beech Be30 
1031 623 N302NC Twin propellers 
King Air 
Single 
1047 1581 Cessna 182 N826DD 
propeller 
1050 1796 Learjet J45 WDR 989 Twin jet 
1052 1925 Beech Be9L N3262R Twin propellers 
Embraer 
1055 2120 SKW5558 Twin propellers 
EM120 
Cessna 750 
1059 2300 EJ947 Twin jet 
Citation 
1102 2480 Beech Be 40 OPT 417 Twin jet 
Embraer XJ 
1106 2780 BTA98 Twin jet 
145 
1118 3427 Beech Be9L N6077X Twin propellers 
Table 2. | Summary of Recorded Target Events 
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B. ANALYSIS OF BACKGROUND NOISE 
1. Power Spectrum of Background Noise 


The background noise level was determined by processing a 300-second 
recording of the ambient environment noise in the absence of any airplane 
passage through the~ relevant functions in MATLAB _ program 
‘auto_Doppler_detector.m’. It is observed that the background noise level is 
highest in the lower frequencies and falls off rapidly as frequency increases. 
Contributors to this low-frequency noise include, but are not limited to equipment 
vibration, wind noise, noise from the ocean to the north, vehicles traveling on 
Highway CA 1 south of NPS, birds, and electronic noise. From the zoom in plot 
in Figure 14, it is seen that the noise spectrum level decreases most rapidly (-20 
dB ) in the frequency range between 0 and 200 Hz. 
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Figure 14. Power Spectrum of Background Noise 
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2. Threshold Determination 


The background noise level is computed to have a mean value of 
=0.0229 and standard deviation of o=0.0058. These values are expressed in 
arbitrary units of pressure squared. This mean and standard deviation are used 
as the baseline for the initial detection threshold. To determine an appropriate 
cross-correlation, R,y threshold, synthetically generated airplane “events” with 
various velocities are cross-correlated with the background noise at various time 
periods and the results averaged. This is repeated for different frequency bands. 
The mean and standard deviation of the peak of the cross-correlation with noise 
are designated as py, and o,. respectively. 


C. ANALYSIS OF DATA — PROPELLER AIRPLANES 


Doppler shifts in the 50-1500 Hz region as the airplane passed through 
CPA are clearly seen for all six propeller airplanes. In most cases, 
measurements of the upper and lower frequencies of the tonals with good SNR 
give a good estimate of the velocity through the use of Equation 2.12. Similarly, 
the rate of change of frequency at CPA is also measured. Processing of the 
tonals is restricted to those having good SNR at the asymptotes before and after 
CPA. In one of the cases, there is an observed change in engine revolutions per 
minute (RPM) just prior CPA, identified by a “knee” in the LOFARgram. The 


velocity in this case was estimated through measurement of fo and frower instead. 


Program ‘auto_Doppler_detector.m’ is run with time domain threshold 
settings of u+150 and Ry, threshold settings of u+0.50;. A total of 131 peaks 
(energy exceeding threshold) and 10 airplane contacts with Doppler shifts are 
detected. A summary of the results obtained with other threshold settings is 
presented in Section E of this chapter. 
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1. Event 1 — N6800M at Time 1021 (74s) 

N6800M is a four-seat, single-engine CESSNA 182P SKYLANE. Used 
mainly for training and general aviation, the C182 has a cruise speed of 150 
Knots Indicated Airspeed (KIAS) or 75 m/s and a stall speed of 49 KIAS (24.5 
m/s). Typical speed during approach for landing is 70 to 80 KIAS (35 to 40 m/s). 









































Tonal 1 Tonal 2 Tonal 3 Program 
to 75s 74s 
fiower | Ffupper | flower | fupper | Flower | Ff upper -- 
67.14 | 91.55 | 97.66 | 134.3 | 195.3 | 274.7 -- 
Af 12.21 18.32 39.70 -- 
fo 79.35 115.98 235.00 122.07 
df/dt -4.91 -6.46 -15.53 -7.761 
v (m/s) 53 54 58 48 




















Table 3. Observed Parameters of N6800M 


There is good correlation between the results obtained by the program 


and manual analysis. 
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Plots of Contact Detected at 74.71 s - N6800M LOFARgram 


Pressure vs Time Plot 
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Figure 15. Plot of Results from Analysis of CESSNA 182 SKYLANE, N6800M 
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2. Event 2—N302NC at Time 1031 (623s) 


The second propeller airplane, N302NC, is a twin-engine, HAWKER 
BEECH SUPER KING AIR 300 BE30/G. A multi-purpose aircraft, the BE 30 is 
used for both civil and military purposes. With a cruise speed of 290 KIAS (145 


m/s) and a stall speed of 75 KIAS (37.5 m/s), the typical approach speeds is 100 
to 120 KIAS (50 to 60 m/s). Only one of the tonals (~ 70Hz) is sharp enough for 


use in the estimation of Doppler shift. 


The other tonals at higher frequencies 


may be “blurred” by both the engines operating at slightly different RPM settings. 
































Tonal 1 Program 
to 623 s 623 s 
f tower f upper - 
85.45 122.11 -- 
Af 18.33 -- 
fo 103.78 94.60 
df/dt -6.35 -4.6566 
v (m/s) 61 53 
Table 4. | Observed Parameters of N302NC 


The results obtained by the program agree with those obtained through 


manual analysis. 
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Pressure vs Time Plot 
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Figure 16. Plot of Results from Analysis of BEECH SUPER KING AIR 300, BE30/G, N302NC 
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3. Event 3 — N826DD at Time 1047 (1581s) 
N826DD is a CESSNA 182T SKYLANE single-engine light airplane similar 


in performance to N6800M. Multiple tonals with Doppler shifts are observed. 
Estimation of velocity is obtained with three tonals. 






























































Tonal 1 Tonal 2 Tonal 3 Program 
to 1581s 1582 s 
fiower | Fupper | fiower | fupper | fiower | f upper -- 
91.55 | 128.2 | 128.2 | 189.2 / 213.6 | 311.3 -- 
Af 18.33 30.50 48.85 -- 
fo 109.875 158.70 262.45 137.33 
df/dt -6.78 -10.85 -16.37 -9.313 
v (m/s) 57 66 64 61 
Table 5. Observed Parameters of N826DD 


A fourth tonal with fo = 855.6Hz is visible. The velocity is estimated to be 
61 m/s based on this tonal. 


The time of CPA is identified by the program to be later than observed. 
This is due to the fact that low-frequency noise (marked with an “X” on the 
LOFARgram) appeared after CPA that skews the total overall energy level. 


The estimate of CPA range for this target is 381 m. This is consistent with 


observations as the airplane was on visual approach to runway 10L. 
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Figure 17. 
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Plot of Results from Analysis of CESSNA 182, N826DD 


4. Event 4 — N3262R at Time 1052 (1925s) 


N3262R is a twin-engine BEECH KING AIR BE/90. Slightly smaller than 
the Super King Air 300 in Event 3, the BE/90 is also used extensively by both the 
military and in civil aviation. With a cruise speed of 260 KIAS (130 m/s) and a 
stall speed of 78 KIAS (39 m/s), typical approach speed is 100 - 120 KIAS (50 to 
60 m/s). Three tonals are observed with only one meeting the required criterion 


for a good estimation of velocity. 
































Tonal 1 Program 

to 1925s 1924s 

f lower f upper _ 

103.80 158.70 -- 

Af 27.45 

fo 131.25 137.33 
df/dt -16.31 -10.70 

v (m/s) 72 48 





Table 6. 


Observed Parameters of N3262R 


From the Ryy vs. Velocity plot, the velocity 48 m/s has the highest peak 
and a second velocity peak exist at 77 m/s. This could be due to interference 
from the noise in the vicinity of f =1930s and a ‘line” near 120 Hz between t = 
1915s to 1920s (Marked with “X” on plot). 


40 





Time (s) 


Fismaupes Tine Cut Plot of Contact at 1924.68s - N3262R LOFARgram 


Pressure 














10 
P 22 oo: 1930 
© i905 
0 ecesecece fe 
= 
B poseeeeneee : seseeeesses =: see ewsaeees Sisvdeetere : wecceiaeeel Beseeseseees ee ee aeveeetees 1920 
1900 1905 1910 1915 1920 1925 1930 1935 1940 1945 0 200 400 600 800 1000 1200 
Time (s) Frequency (Hz) 
Tonal Analysed Plot of Rxy vs Velocity 
32 T 
Threshold : 
30 1 . 
1930 
28 ooo ccemecccneecerees 
1925 Ee pe RE sn Pc sienn canes vane et Wiesuts toms eces Bence oranare dosebin tovniacce Neate etvteeseases 
24 af ccccccccccnceecccncccnsecsagesscssccsa ess g esses Faneassesgessssssssessssgesesssssssssmbassessssssssses 
1920 
22 Sacer eeeeeenaeeeeseees nese egesseeseesee esse esseeeseess nes gesesssessessnsaesssssssssesssmssessbasssees 
1915 20 : : H : : : 
100 120 140 160 180 30 40 50 60 70 80 90 100 


Frequency (Hz) Velocity (m/s) 


Figure 18. Plot of Results from Analysis of BEECH KING AIR BE90, N3262R 
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5. Event 5 — SKY5558 at Time 1055 (2120s) 


SKYWEST Flight 5558 departed Monterey for San Francisco, flying over 
the tower at Herrmann Hall at a higher altitude than the other airplanes. The 
airplane is a twin-engine turboprop Embraer EMB-120 Brasilia. A commuter 
airplane designed for commercial use, it has a passenger capacity of 30, a cruise 
speed of 300 KIAS (150 m/s) and a stall speed of 110 KIAS (55 m/s). The 
airplane was on her climb out during the recording. One tonal line that was 


clearly distinct and another “half tonal” centred on 240.9 Hz were both used for 



































estimation. 
Tonal 1 Tonal 2 Program 
to 2121s 2119s 
f tower | fF upper Tower Tupper és 
100.7 | 143.4 201.1 - a 
Af 21.35 39.8 -- 
fo 122.05 240.9 131.23 
df/dt -6.1833 -10.58 -6.208 
v (m/s) 60 57 50 














Table 7. Observed Parameters of SKY5558 


The results are in agreement with the ground truth. Estimated CPA time is 
earlier than observed, due in part to the high level of low frequency noise at time 
2119s (marked on plot). 
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Figure 19. Plot of Results from Analysis of EMBRAER EM120, SKY5558 
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6. 


Event 6 — N6077X at Time 1118 (3427s) 


The last of the propeller airplanes, N6077X is also a twin-engine BEECH 
KING AIR 90. There is an observed increase in the engine RPM just before the 


airplane passes through CPA, which is seen in the shift (knee) in the LOFAR 


lines. 


airplane’s estimated parameters are as tabulated. 


CPA time is estimated from the sound pressure — time plot and the 





















































Tonal 1 Tonal 2 Tonal 3 Program 
to 3425 s 3425 s 
f tower | Fupper | f tower fupper | Fiower | f upper -- 
140.4 - 152.6 - 225.8 - -- 
Af 24.40 36.60 55.00 -- 
fo 140.40 189.20 280.80 283.81 
df/dt -5.43 -5.78 -16.27 -6.208 
v (m/s) 60 66 67 30 
Table 8. | Observed Parameters of N6077X 


Although the passage is correctly identified by the program, the R,y value 


is barely above threshold. 


The change just prior to CPA and the resulting 


absence of a complete Doppler shift gives a poor estimate of df/ot and velocity. 
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Figure 20. 


Plot of Contact at 3424.68s - N6077X 


LOFARgram 
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Plot of Results from Analysis of BEECH KING AIR BE 9L, N6077X 


D. ANALYSIS OF DATA — JET AIRPLANES 


The LOFARgrams of the jet airplanes are observed to be broadband with 
energy extending into the 12 kiloHertz (KHz) region. While the targets are 
detected as “peaks” in the time domain, there are however no distinct tonals that 
can be tracked. The Doppler shift patterns that are visibly observable for the 


propeller airplanes are absent here. 


A total of four jet airplanes are recorded. All of them have twin engines 
and cruise speeds of 450 to 600 KIAS (225 to 300 m/s). Approach speed is 
typically 150 to 200 KIAS (75 to 100 m/s) for these airplanes, which is 
considerably higher than the propeller airplanes. 


46 


1. Event 7: LEARJET J45 (WDR989) at Time 1050 (1796s) 


Sound Pressure vs time plot - Learjet 45, WDR989 
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LOFARgram (Learjet 45, WDR989) 
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Figure 21. Pressure vs. Time Plot and LOFARgram of LEARJET J45, WDR989 
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2. Event 8: CESSNA Citation C750 (EJ 947) at Time 1050 (1796s) 





Sound Pressure vs time plot - Cessna C750, EJ947 
8 1 
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Figure 22. Pressure vs. Time Plot and LOFARgram of CESSNA C750, EJ947 
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3. Event 9: BEECH Be40 (OPT417) at Time 1102 (2480s) 





Sound Pressure vs time plot - BE40, OPT417 
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Figure 23. Pressure vs. Time Plot and LOFARgram of BEECH Be40, OPT417 
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4. Event 10: EMBRAER XJ 145 (BTA98) at Time 1106 (3427s) 


Sound Pressure vs time plot - Embraer ERJ 145, BTA98 
3 T T T T ij 
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Figure 24. Pressure vs. Time Plot and LOFARgram of EMBRAER ERJ 145, 
BTA98 
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E. PERFORMANCE OF THE AUTO-DETECTION PROGRAM 


Using the threshold settings u+150 and R,y threshold setting of /-+0.50, , 
the algorithm successfully identified all six propeller airplanes. Of the 10 
“detected with Doppler calls,” there were four instances of false alarm. 


The processing is repeated with other combinations of threshold settings 
of u+To and Ry, threshold setting of u-+Ro;, The results are tabulated below. 
The output analysis reports are included in Appendix H. It is seen that the 
energy threshold setting affects mainly the number of peaks detected while the 
Ry, threshold setting affects mainly the number of Doppler shifts observed. 
Above a certain threshold setting, however, the propeller planes are “cut off’ and 


the number of correct detections falls off. 




















Threshold Ryy threshold No. of No. of 
Run Remarks 
Setting, T setting, R Peak Doppler 
1 15 0.5 131 10 4 false contacts 
N6077X not 
2 15 1.0 131 5 
detected 
3 20 0.5 113 7 2 false contacts 
N6077X not 
4 20 1.0 113 5 
detected 
Only 4 
airplanes 
5 48 1.0 65 5(4) 
detected. 1 
double-count 




















Table 9. Summary of Results using auto_Doppler_detector.m 
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It is clear from the results that the correct identification of the CPA time fo 
impacts the results significantly. In Events 1 and 2, to determined by the 
algorithm is close to the observed CPA time and the estimates of velocity and 
range are close to the ground truth. In addition, the Ryy is significantly larger than 
the threshold value and the R,y -velocity curve has a sharp peak with steep 
sides. In Event 3, we see the “rounding” off of the peak as CPA time starts to 


deviate from observed. 


As the experiment was conducted using targets of opportunity rather than 
controlled contacts with known velocity and passing distance, an accurate 
assessment of the algorithm’s performance and accuracy cannot be made. 
However, the velocities obtained are all within typical velocities for the respective 


types of airplanes on approach to landing. 


As observed in the various LOFARgrams of jet airplanes, the signals are 
broadband and spread over a wider range of frequencies than the propeller 
planes. The absence of distinct tonals enables the algorithm to correctly 
discriminate between these targets and the propeller planes using the Doppler 
filter. 


F. CONCLUSION 


The experimental results show that using the presence or absence of the 
Doppler shift enables selection of a low threshold for energy detection without 
concern that the detector will be overwhelmed by false alarms. The number of 
initial detections using energy criterion is filtered down to a manageable number 


by the secondary criterion involving Doppler. 


It is seen that the Doppler shift pattern is useful in two main ways. First, 
even when the threshold for energy detection is set at a relatively low level, the 
use of a “second level filter” significantly reduces the number of contacts that 
have to be analysed by automatically filtering out those that do not meet the 


defined criteria. This reduces P;, when a false alarm is taken to be a 
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misidentification as well as a “mis-detection.” Finally, valuable information 
(velocity and CPA range) is also derived automatically although these estimates 


may be limited in accuracy. 


53 


THIS PAGE INTENTIONALLY LEFT BLANK 


54 


V. CONCLUSION AND RECOMMENDATIONS FOR FUTURE 
DEVELOPMENTS 


A. ACCURACY OF RESULTS 


The use of the Doppler shift in the decision process has been shown to 
improve the performance of an automatic detector for a target passing through 
CPA compared to detection based solely on energy levels. Doppler shift as an 
additional criterion allows the threshold to be set near the noise floor while 
achieving a lower Pr. The performance of the detector in the estimation of 


contact velocity and range relies on the following factors. 


1. Determination of CPA Time 


The time of CPA is the basis for the analysis of the Doppler pattern. In our 
study, the CPA time is based on the assumption that the total energy in the 
frequency domain peaks during CPA. This is subject to two possible ways in 


which errors affecting the results can be introduced. 


First, the Doppler-Range (frequency-time) ambiguity due to the Fourier 
transform process places a limitation on the time resolution that can be obtained 
for a specific desired frequency resolution. The CPA time that is determined 
within the frequency domain is hence an indication that CPA is happening 
sometime within a time block rather than at a specific resolved time. Depending 
on the speed and CPA range, this could translate to significant distance between 
the measured and actual CPA. Secondly, the possibility of interfering contacts or 
background noise occurring in the vicinity of CPA time can displace the peak and 


hence affect the accuracy of the results. 


While the latter is due to factors that are beyond control, improvements 
can be made to mitigate the effect of the former. Determining the CPA time in 
the time domain by tracking the trend of the magnitude of the received signal and 
tagging CPA to the maximum will allow CPA time to be declared with a better 
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resolution. Alternatively, since it is the total energy across all frequencies, the 
power spectral density could be calculated with a short integration time, 
sacrificing frequency resolution, purely to estimate the CPA time. Finally, if the 
sensor has other information (e.g., contact bearings), these could be used to 
refine the CPA time. 


2. Measuring the Rate of Frequency Change, df/dt 


The rate of frequency change, df/dt is one of the two parameters that 
generates the template Doppler curves used in the cross-correlation with the 
measured signal. Once determined, this is not varied. Only the velocity is varied 
in the cross-correlation process. The df/dt is estimated based on the assumption 
that in the time blocks before and after CPA, the strongest frequency 
corresponds to that of fo. This again could be subject to errors from interfering 
noise sources. In addition, the number of differences that could be used in 
determining df/dt is constrained by the need to limit the analysis to the vicinity of 
CPA to ensure that the assumption in Eq. 2.8 holds. This, in turn, is determined 
by the passing range and the velocity of the contact. A fast contact at close 
range means that the contact is within the acceptable range of values near CPA 
for only a very short period of time and the number of points used to determine 
df/dt is small. 


An improvement might involve the use of cross-correlation to determine 
df/dt instead. Lines of different slopes could be correlated to the measured 
signal for the best estimation of df/at. 


3. Number of Tonals 


Our analysis only involves the use of a single tonal even in cases where 
there are several tonals that exhibited good SNR for analysis. A more accurate 
estimate could be obtained if more tonals are used. The results from reading 


several df/dt values could be plotted against their corresponding fp values and the 
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slope obtained with linear regression could be used to provide a better estimate 


y 
of —. 


oO 


B. CONCLUSION 


The objectives set forth for this thesis have been met. It has been 
demonstrated with one hour of recording that the Doppler effect can be used to 
improve the passive automated detection of contacts. The estimates of velocity 
and range for the airplanes are within expectations based on ground-truth 
observations. Finally, the detection algorithm developed here is compact and 
could realistically be deployed in acoustic sensors as part of an underwater 


surveillance network. 


C. RECOMMENDATIONS FOR FUTURE DEVELOPMENT 


It is recommended that steps to improve the determination of CPA time 
and df/dt be explored in the future. Improvements in the estimate of these 
parameters will improve the overall detection performance in terms of Py and Pra. 
In addition, the current detection algorithm could be extended to use more than 


one tonal. 


Another follow-up step is to extend the detection algorithm to the 
determination of contacts in an underwater environment. This extension will be 
complicated by the complexity of underwater acoustics, including multipath 
propagation, refraction, and noise variability. However, underwater detection will 


be aided by the slower target speeds and longer integration times. 
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APPENDIX A. TECHNICAL SPECIFICATIONS 


Technical Specifications for Data Acquisition Module 


USB Dynamic Signal Acquisition 





NI USB-9233 


© 24-bit resolution Recommended Software 
© 102 dB dynamic range © LabVIEW 

© 50 kS/s max rate per channel © LabVIEW Sound and 

© 4 simultaneous anglog inputs Vibration Toolkit 


@ +5V input range . 

© AC coupled with IEPE power © LabWindows/CVI 

© Hi-Speed USB 2.0 © Measurement Studio 
© DAQmx 


LabVIEW Order Analysis Toolkit 








Overview 


The National Instruments USB-9233 is a four-channel dynamic signal 
acquisition module for making high-accuracy measurements from IEPE 
sensors. The NI USB-9233 delivers 102 dB of dynamic range and 
incorporates IEPE signal conditioning for accelerometers and 
microphones. The four USB-9233 input channels simultaneously acquire 
at rates from 2 to 50 kHz, In addition, the module includes built-in 
antialiasing filters that automatically adjust to your sampling rate. 


Hardware 


Data rates on the four simultaneous input channels of the USB-9233 range 
from 2 to 50 kHz. Each signal is AC coupled, buffered, analog prefiltered, 
and sampled by a 24-bit delta-sigma analog-to-digital converter (ADC) that 
performs digital filtering with a cutoff frequency that automatically adjusts 
to your data rate. The USB-9233 module features a voltage range of +5 V 
and a dynamic range of over 100 dB. Communication to the USB-9233 is 
Hi-Speed USB 2.0, guaranteeing data throughput. 

The USB-9233 is a nonisolated module that provides +30 V of 
overvoltage protection (with respect to chassis ground) for IEPE sensor 
connections. IEPE excitation and AC coupling are not software-selectable 
and are always enabled. 

The USB-9233 uses a method of A/D conversion known as delta- 
sigma modulation. If the data rate is 25 kS/s, each ADC actually samples 
its input signal at 3.2 MS/s (128 times the data rate} and produces 
samples that are applied to a digital filter, This filter then expands the 
data to 24 bits, rejects signal components greater than 12.5 kHz (the 
Nyquist frequency}, and then digitally resamples the data at the chosen 
data rate of 25 kS/s. This combination of analog and digital filtering 
provides an accurate representation of desirable signals while rejecting 
out-of-band signals. The built-in filters automatically adjust themselves 
to discriminate between signals based on the frequency range, or 
bandwidth, of the signal. 


Software 


The USB-9233 uses NI-DAOm«x as the hardware and OS interface, so you 

can build automated test systems or integrate the USB-9233 with other 

hardware, including National Instruments multifunction DAO products 

and CAN devices, through NI-DAQm« function calls for NI LabVIEW, 

LabWindows/CVI, and Measurement Studio. 

Other useful NI-DAQ driver software features include: 

© Automatic Code Generation — DAQ Assistant is an interactive guide 
that steps you through configuring, testing, and programming 
measurement tasks and automatically generates the necessary code 
for LabVIEW, LabWindows/CVI, or Measurement Studio. 

¢ Test Panels — With NI-DAQ, you can test all of your device 
functionality before you begin application development. 

© LabVIEW Integration — All NI-DAO functions use the waveform data 
type, which carries acquired data and timing information directly into 
more than 400 LabVIEW built-in analysis routines for display of results 
in engineering units on a graph. 


Analysis Software 

Because of its bandwidth, the USB-9233 is well-suited for applications 
in noise and vibration analysis. These applications are specifically 
addressed in the LabVIEW Sound and Vibration and LabVIEW Order 
Analysis toolkits. By using these toolkits in conjunction with a USB-9233 
module, you can produce power spectra, frequency responses, fractional- 
clave analysis, sound-level measurements, order spectra, order maps, 
order extraction, and sensor calibration. The USB-9233 is ideal for a wide 
variety of mobile/portable applications such as industrial machine 
condition monitoring and in-vehicle noise, vibration, and harshness testing 


“Winstruments 


(From: National Instruments, n.d.) 
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USB Dynamic Signal Acquisition 












































































Sp ecifications Dynamic characteristics 

The following specifications are typical for the range 0 to 60 °C unless Passband 

otherwise noted. ‘s aud ceuantal sa 

ss $25 kS/s 045 - fs 0.06 B 045 - fs 334° 

Input Characteristics SIKSIs 042. fs 006 Oat. ts 2 

Number of channels... 4 analog input channels 

ADC resolution.......... 24 bits . a Stophand aprsad Alizs-Free 

TYPO OF ADC ssc ccccsccsssecscseceriessccessssins Delta-sigma oa <25 KS's TR _— ont — - 

(with analog prefiltering} >I5KS/s 068 92 4B 64 fs 032. fs 

Data rate {fs} ——— —— —————s 
Minimum 2 kS/s 
Maximum .. 50 kS/s Crosstalk 

Master timebase {internal} Paired channels 
=a bs on ‘nahin (0 and 1, 2 and 3} 100 dB at 1 kHz 
reek i 100 7 Nonpaired channel -110 dB at 1 kHz 

i aed a Den mes Commen-mode rejection ratio (CMRR} 

A cacti 9 Minimum... A dB, fin <1 kHz 

Oth Arequency ; Kpical....... 80 dB, fin< 1 kHz 
: ye : \yp Spurigus-free 
Wi eheee tke a rE ME dynamic range (SFR)... 120 dB (fin = 1 kHz, -60 dB FS} 
Voltage Tull-scale range ; Idle channel noise and ngise density 
Typical ....... 5.4 Voeak ‘ais ame Hise ss 
Mini 5Y, Channel 
pikes 5 BV. Noise & BFS $B dB FS 102 dB FS 
7 eS Noise density 400 nvAVHe 400 nvAVH2 900 nVANHE 

Common-mode voltage 
(Al- to earth QFOUNG) 0. seessesessseoes +2¥ Input impedance 

IEPE excitation current Differential (AC)... $300 kQ 
Minimum... 2.0mA Al- (shield) to chassis ground 502 
Typical ....... 22mA Distorti 

; istortion 

IFPE compliance voltage 19¥ max Harmonie {THD} 

Overvoltage protection (with respect to chassis ground} 1 id, | 10 kHz, 10kHe, 
For an IEPE sensor connected lee een e BRS me Sahat 
to Al+ and Pile scssctesdesieicnnnnnnesnsnct 230V 20 dB FS 95 dB 9048 oe 
For a low-impedance source connected ea a ee ee 
to Alt and Ale eee =O tO 30 V : : 

Input delay {in seconds} mares a oe me 
een ie 4:1 amplitude ratio... 80 dB 
ey ee et CCIF 11 kHz/12 kHz 

Accuracy (0 te 60 °C) 1:1 amplitude ratio... -93 dB 

Error Accuracy 
Calibrated max +0308 
Calibrated typ 40.1 6B 
Uncaliteated max 106 dB 

Accuracy drift 
Typical ....... 0.001 dB/°C 
Maximum ...... 0.0045 dBc 

Channel-to-channel matching 
Gain 

Maximum... . 027 dB 
Typical....... 0.07 dB 
Phase (fin in kHz}... fin 0.077° + 0.067° 





BUY ONLINE at ni.com or CALL (800) 813 3693 (U.S.) 


(From: National Instruments, n.d.) 
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B. Technical specifications of ACO Pacific Microphone 


1. Cartridge model 7046 
ACO JAPAN ACO Jt 
CONDENSER MICROPHONE E — - 
FREQUENCY RESPONSE r+ 
CARTRIDGE TPE 7046 "9 es ER 
SERIAL No, 36930 





“253 Bre.) V ver pascal or 54.33 mV ver vance *'” | 


Open Circuit Correction Factor: 

Ke= a 
Cartridee Capadtance: 

c= 205 oF 


Lockace Resistance testes at S2%relative humidity 





























2. 4012 — 1/2 Inch Preamplifier (From: Aco Pacific, n.d.) 
Specifications 


. 28 Vdc, 22 pF (DM2-22 Dummy Mic) Terminated 
. Frequency Response: 

<2 Hz to >200 kHz +/- 0.5 dB, 

0.5 Hz to >> 200 kHz +/- 2 dB 
. Insertion Loss: <<0.25 dB @ 1 kHz 
. Noise: 

"A" Weighted - <1.2 uV typ. 

Linear 20 to 20 kHz - 3 uV 


. Input Impedance: >50 G//0.1 pF 
. Output Voltage: >7 Vrms @28 Vdc Typical 
>15 Vrms @50 Vdc 
. Power Requirements: 
28 Vdc @ <1 mA (HP option <2 mA) 
50 Vdc @ <2 mA (HP option <4 mA) 
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APPENDIX B. MATLAB CODE : AUTO_DOPPLER_DETECTOR.M 


% Auto Doppler Detector Program. 

% Purpose : This program search a time domain input signal for Doppler 

% Shifts and computes the velocity and passing range of the contact. 

% Detection is based on user determined threshold. 

% Input: O07 _aug_08.txt : Input time domain data file 

% threshold.mat : threshold data generated using 

% ‘generate_threshold.m' 

% Output: 07 _aug_08_result.txt: Tabulated time, frequency, velocity % 
and range results from analysis. 

% Various plots of detected peaks with Doppler Shift. 

% Variables : 

% 

%_ infile, outfile : File names of input and output file. 

% Fs : Sampling frequency derieved from input data. 

% time_length (s) : User specified length of data to process in each 

% cycle. Recommended to note MATLAB limitations 

% of array size. 

% time_duration(s) : User specified expected total duration of each 

% passage. 

% input : Time domain signal data read from infile. 

% Y : input data tagged with time of data information. 

% del_f,int_t : frequency resolution desired and integration 

% time for computation of power spectral. 

% threshold_setting,Rxy_threshold : User specified threshold. 

% Functions called: time_to_freq.m 

% find_peak.m, 

% calc_df_dt.m 

% peak_compare.m, 

% generate_template.m (used by function peak_compare.m) 

% 

% Comments : Input text file is generated using NATIONAL LAB LABVIEW 

% Express software, converted from .tdms to .txt. The 

% first 7 lines contains setting information related to 

% the recording and data starts on the 8th line. 

% : Desired threshold is through the setting of variable 

% "threshold_setting". This adds the mean threshold the 

% corresponding number of 'Standard Deviation’. 

% : Output file name is tagged to input file name. the suffix '_results' is added. 


% Written by How Whye Keong, NPS, 05 Jan 2009 


WT % V0 % Yo %o %o Vo Vo YoYo %o Yo Vo %o Vo Vo Vo Yo Yo Yo Yo Yo Vo Vo Vo Yo Yo Yo Yo %o Yo Yo Vo Vo Vo Yo Yo Yo Yo Yo Yo Vo Vo Vo 
63 


close all 
clear all; 


Vh% Setting of Variables %%% 


infile='07_aug_08.txt’; % Name of input file 

time_length = 300; % length of each segment to be read in sec. 
time_duration = 20; % length of the expected passage of interest 
del_f = 4; % Desired freq resolution. 

int_t = 1; % \ntegration time 


threshold_setting = 20; 
Rxy_threshold = 0.5; 


© Desired thershold setting. 
Desired Rxy threshold setting. 


se 3k 


fid=fopen(infile,'r’); % Open file for reading only 

line1 = fgetl(fid); % Reads line 1 as character string 
line2 = fgetl(fid); % Reads line 2 as character string 
line3 = fgetl(fid); % Reads line 3 as character string 
line4 = fgetl(fid); % Reads line 4 as character string 
lineS = fgetl(fid); % Reads line 5 as character string 
line6 = fgetl(fid); % Reads line 6 as character string 
line? = fgetl(fid); % Reads line 7 as character string 


start_time=line4(1:20); 

delta_time=line6(1:8); 

delta_time = str2double(delta_time); 

Fs = round(1/delta_time); % Sampling Frequency. Round to ensure integer 


clear line1 line2 line3 line4 lined line6 line7; 


% Read the stipulated data from input file and tag the corresponding 
% real time to data points. 


% The output is matrix y of sound pressure of time_length with padding 
% of buffer_length to the begining and end. Exception is the first 
% 


o segment. 
time_tracker = 0; % Initialise time_tracker, peak_count & 
peak_count = 0; % detection_count. 


detection_count=0; 
time_duration = floor(time_duration/2); 
for j = 1: 100 


DhW%% Read data from file %%hhhHhV%%%%%%%%o%% Vo Vo % Vo %o Vo %o% Vo %o%o Vo Vo %o Vo %o 
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% Determine number of data points to read 


if j == 
data_pts= (time_length + time_duration) * Fs; 
else 
data_pts= time_length * Fs; 
end 
disp("’); 
disp(‘Loading data ...'); 
input = fscanf(fid,'%f\n',[1 ,data_pts])’; % Reads data from file 


if isempty(input) 
disp('End of file reached.'); 
break; 

else 
% Create real time vector (Required for plotting) 
time_vec = linspace(1, length(input),length(input))’; 
time_vec = (time_vec - 1) * delta_time; 
time_vec = time_vec + time_tracker; 


time_tracker = (j* time_length) + time_duration; 
y =[time_vec input]; 


% Store the last segment of data to be attached to the next 
% data segment read. 
if j == 
y_buffer_before = y(((time_length - time_duration)*Fs)+1 : ... 
length(input) , :); 
else 
y = [y_buffer_before ; y]; 


m = size(y); 
y_buffer_before = y ((m(1)-(2*time_duration*Fs))+1 : m(1), :); 
end 
end 


disp(['Data from segment ',num2str(j),’ loaded...']); 


%*%% Convert to Frequency domain and locate energy peaks %%% 
% Performs FFT on input signal 

disp (‘Performing FFT...'); 

[f,t,lofar] = time_to_freq(y(:,2),Fs,del_f,int_t); 


t=t+y(1,1); % Align time vector of LOFARGRAM to real time 
disp (‘FFT completed’); 
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%%% Locate the instances where energy is above threshold %%% 
% Load threshold data. 


if j== 
load threshold.mat; 
threshold = threshold + threshold_setting*SD; % Set threshold 
threshold_vec(4,:) = threshold_vec(2,:) + Rxy_threshold * ... 
threshold_vec(3,:); 
peak = zeros(1,9); % Initialise a matrix to store results 
end 


% Locate energy peaks 
disp('Finding peaks...'); 
[p_out] = find_peak(lofar,f,t,threshold,time_duration); 


*V%% Analysis of Peaks %%% 


% Compute the df_dt for each of the peak and use this information 
% to synthetically generated Power Spectra for different 
% velocities. Find best match using correlation. 


[m,n]= size(p_out); 
peak_count=peak_count + m; 


ifm>0 
disp([num2str(m),' peaks found.']); 
disp(’ ’); 


fork=1:m 


disp(['Analysing peak ',;num2str(k),’.']); 

[p_out(k,6),p_out(k,7)]=calc_df_dt(lofar, f, t, p_out(k,:)); 

[p_out(k,8),p out(k,9),p_out(k,7)] = peak_compare ... 
(lofar, f, t, p_out(k,:),time_duration,threshold_vec); 


if p_out(k,7) == 2 
% Plot Pressure-time graph 


detection_count = detection_count + 1; 
start_time_I = find (y(:,1) >=... 

(p_out(k,2) - 2*time_duration),1,'first'); 
end_time_l = find (y(:,1)> ... 

(p_out(k,2) + 2*time_duration),1,'first’); 


subplot(2,2,1); 

plot (y (start_time_I :end_time_I ,1).... 
y(start_time_l:end_time_I,2)); 

xlabel('Time (s)’); 

ylabel(‘Sound Pressure (AU)'); 

grid on; 
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title (‘Sound Pressure vs Time Plot’); 


clear start_time_| end_time_|; 
end 
end 
end 





% Copy result to Master Table 


if j== 
peak = p_out; 
else 
peak = [peak ; p_ out]; 
end 
end 
fclose(fid); % close ASCII file 


% Generate output report 


flag_interpret = char('No Doppler Shift. df/dt inadmissible... 
‘No Doppler Shift Detected. Rxy less than Threshold... 
‘Contact with Doppler Shift found.'); 


outfile=[infile(1:9),'_results.txt’]; % Name of output file 
fid=fopen(outtile,'w'); 


forintf(fid,'Results from analysis of %s.\n\n',infile); 
fprintf(fid," Threshold settings : %6.2f.\t Rxy_threshold : %6.2f.\n', 
,Rxy_threshold); 


fprintf(fid,"Threshold for energy based detection = ... 
%7.Af.\n' threshold); 
forintf(fid,'No of peaks found = %4.0f.\n',peak_count); 
fprintf(fid,'No of Doppler Shifts detected = ... 
%3.0f.\n\n',detection_count); 

forintf(fid,'Start time : %s.\n\n',start_time); 
fprintf(fid,"Peak No. Time Freq Result 

Velocity Rxy Range \n‘); 


for i = 1: peak_count 


fprintf(fid,"%3.0fMt \t %7.2At %7.2ft %s\t '.... 
peak(i,2),peak(i,3),flag_interpret((peak(i,7)+1),:)); 


if peak(i,7)== 
R_0 = -((peak(i,8)*2)/(peak(i,6)))*(peak(i,3)/343); 


fprintf(fid,'"%7.2ft %7.2f\t %7.2F \n',peak(i,8),peak(i,9),R_O); 
else 
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threshold_setting 


forintf(fid,"\n'); 
end 


end 
fclose(fid); 
disp([Analysis Completed. Results written to file ",outfile,""]); 


PWM PWWWWWV% End of Program %%%%% %%%%%% Vo % Vo Vo %o%o% Yo Vo %o Yo %o 
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APPENDIX C. MATLAB FUNCTION CODE : time_to_freq.m 


function [f,t,lofar] = time_to_freq(y,Fs,del_f,int_t) 
Computes the Power Spectral Density. 


3 


% 
% 
% 


Purpose : This function performs the Fourier Transform of the 
time-domain input y and returns the Power Spectral Density estimate 
via Welch's method. 


te > ee * 


% \nput Variables 

% Y_—: Signal in time domain. 

% Fs : Sampling Frequency. 

% del_f : Desired Frequency spacing (freq resolution) 
% int_t : Integration time 


° 


% Output Variables 


% ~ lofar : Matrix of Power Spectrum . 

% fF : Frequency Matrix for ‘lofar' (x-axis) 

% t : Time matrix for ‘lofar’ (y-axis) 

% Local Variables 

% NFFT : Number of FFTs to perform 

% avg : No of averages. 

% chunk : Length of each time domain data segment 
% N : No. of time domain data segments 

% Functions called : none 

% Written by Prof. Daphne Kapolka, NPS 

% Adapted by How Whye Keong, NPS, 10 Feb 2009 


LT % V0 % Vo %o Vo Vo %o Yo %o Yo Vo Vo %o Vo Vo Yo Yo %o Yo Vo Vo Vo Vo Yo Yo Yo Yo %o Yo Yo Vo Vo Vo Yo Yo Yo Yo Yo Vo Vo Vo Yo Yoo 


NFFT = 24(ceil(log2(Fs/del_f))); % No. of FFTs 
window = hamming(NFFT); 


time_FFT = NFFT/Fs; 
avg = floor(int_t/time_FFT); % No of averages 
if avg == 
avg = 1; 
end 
chunk = NFFT*avg; %_ length of time data needed to get averages 


N = floor(length(y)/chunk); —=% Total number of chunks available 


lofar = ones(N,NFFT/2+1); % initialize matrix for speed 
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x = ones(1,NFFT/2+1); 
t = (1:N)' * (chunk/Fs); % time corresponding to end of the 
% ith chunk 


fori=1:N % calculate avg power spectrum for each chunk 


low = 1+(i-1)*chunk; 
high = i*chunk; 


[x,f] = pwelch(y(low:high),window,NFFT/2, NFFT,Fs); 


X =x" 


lofar(i,:}=x; |% make the ith row of the lofargram the power 
% spectrum for that chunk 


end 


DWV %%% 0% %%%%o%%%% End of function %%%%%% %%%%%%%o%o% 
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APPENDIX D. MATLAB FUNCTION CODE : find_peak.m 


function [peak] = find_peak(lofar_func,f_func,t_func,threshold ,time_duration) 
% Search for peaks in the frequency domain of input signal. 


% Purpose : This function search for energy peaks in a signal in the 

% input frequency - time matrix (lofar_func). This is done by summing 
% the total energy across all frequenies within a time band and 

% compare it to a user determined threshold. The function also search 
% for the tonal (frequecy) with the highest magnitude above a preset 
% frequency. 


% The function's output is a matrix (peak) containing the time (peak 
% time) when the threshold was exceeded and strongest tonal's 
% frequency. 


% Input Variables 

%  \ofar_func: signal in frequency domain. 

% f_func : frequency matrix of signal corresponding to % 
lofar_func. 

% t_func : time matrix of signal corresponding to lofar_func 

% threshold : The threshold for peak detection. 

%  time_duration : The length of passage. Used to determine how the 
% length of buffer segment. 


% Output Variables 
% peak : Matrix with time(peak) when threshold was exceeded. 
% Structure of Matrix 'peak' 


% PeakNo. Time Frequency Time_index Frequency_Index 


Of | ete ee ete ee ee 
% Local Variables 

% threshold : Threshold signal compared to 

%  freq_min : Lower boundary to search for strongest tonal. 

% energy : Energy contained in a time band 

% num_peak : Count of instances where threshold is exceeded. 
% *.* I : Index of parameter *.*. 

% a,yn,i : dummy variables and loop index. 


% Functions called : none 


% Comments : 

% Since the input time domain signal is padded at both ends with data 
% of length = time_duration, any peaks found in these regions is 

% omitted from the results since these regions overlaps with either 

% the previous or next data segment. 


% Written by How Whye Keong, NPS, 10 Feb 2009 


WT V0 % Vo %o Vo Vo Vo %o%o%o Vo Vo %o Vo Vo Yo Yo %o Yo Vo Vo Vo Vo Yo Yo Yo Yo %o Yo Yo Vo Vo Vo Yo Yo Yo Yo Yo Vo Vo Vo Yo Yo Vo 
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freq_min=90; % This is the lower boundary for the frequency band 
% to search in. 


%%o 

energy = sum(lofar_func,2); % sum total energy over all frequencies 
% ~ figure; 

%  plot(t_func,energy); 

% — title(‘Total Energy vs Time’); 


%  xlabel('Time (s)'); 
%  ylabel('Total Energy’); 


%%o 
% Find start & stop times where signal exceeds threshold and store in 
% matrix 'peak_times' 

over_threshold = energy > threshold; 


% Locate start and stop times 


flag = over_threshold(1 : length(over_threshold)-1); 
flag_2 = over_threshold(2 : length(over_threshold)); 
flag = xor(flag,flag_2); 


peak_times = t_func(flag); 
num_peak = floor(length(peak_times)/2); 
% Tabulate the peak time, strongest tonal Freq and their indexes. 


n_time_seg = floor(time_duration/(t_func(2)-t_func(1))); 
% time index of buffer section 


peak = zeros (1, 5); % Output matrix 
if num_peak > 0 


peak_counter = 1; 
fori = 1 : num_peak 


peak_st_l = find(t_func == peak_times ((2*i)-1)); 
peak_et_| = find(t_func == peak_times (2*i)); 


[a,n] = max (energy (peak_st_I : peak_et_l)); 
% n = index of maxmimum energy in time span 


peak_time_| = peak_st_I+n- 1; 
% \ndex of peak time in t_func 
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% Search for strongest tonal in frequencies from freq_min Hz 


% onwards. 


f_min_| = find (f_func > freq_min , 1 ,'first’); 


[a,n] = max (lofar_func(peak_time_I , f_min_I: ... 


f_max_l=f_min_l+n-1; 





% Only record peaks outside of buffer region 


if (peak_time_| > n_time_seg) && (peak_time_I<... 
length(t_func)- n_time_seg) 





% Peak number 


peak (peak_counter,1) = peak_counter; 


% Time of Peak 


peak (peak_counter,2) = t_func(peak_time_l); 


% Freq of strongest tonal 
peak (peak_counter,3) = f_func(f_max_l); 


% Index 
peak (peak_counter,4) = peak_time_|; 
peak (peak_counter,5) = f_max_]; 


peak_counter = peak_counter+1 ; 
end 
end 
else 
peak = []; 
disp ('No Peaks found’); 
end 


% \ndex of min_freq 


length(f_func))); 


DMM VV VWWPWVV%% end of function %~hhVV%%%%%%%%o%o% 
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APPENDIX E. MATLAB FUNCTION CODE : calc_df_dt.m 


function [df_dt,flag] = calc_df_dt(lofar_func,f_func,t_func,peak_func) 
Computes the rate of change of frequency. 





3s 


% Purpose : This function computes the rate of change of frequency of 
% atonal. "Tracking" of the tonal through time is based on the 

% assumption that within the time period of interest and in the 

% vicinity of the tonal's frequency, the tonal will have the highest 

% magnitude. 


° 


° 


% Input Variables 

% lofar : Matrix of Power Spectrum 

% f func : Frequency matrix of power spectrum 
% tfunct : time matrix of power spectrum 


%  f_func and t_func is used to compute frequency resolution and time 
% spacing data points 


% peak_func : matrix containing information on time and 

% frequency of tonal of interest. 

% Format as follows : 

% 

% PeakNo. Time Frequency Time_index Frequency_Index 
Ofte te wee) 6 hue 9 cen So bakes re ok 


% Output Variables 
% df_dt : The mean rate of change of freq. for tonal of interest 


% Local Variables 


% max_v : Maximum expected velocity. Affects the freq range 
% around tonal in which changes are looked for. 

% del_f, int_t : frequency and time between data points. 

% span : No. of bins based on max_v and del_f. 

% no_point  : No. of time periods before/after tO 

% used to calculate mean df/dt. 


% Functions called : none 
% Written by How Whye Keong, NPS, 11 Feb 2009 


WV V0 % Yo %o Vo Vo %o Yo %o Yo Vo %o Vo Vo Vo Yo Yo Yo Yo Vo Vo Vo Vo Yo Yo Yo Yo %o Yo Yo Vo Vo Vo Yo Yo Yo Yo Yo Vo Vo Vo Yo Yo Vo 


max_v = 100; % Max velocity of contacts expected (m/s) 
c_air = 343; % Speed of sound 
no_point = 1; % Number of points to base df/dt calculations 


% Calculate freq and time spacing. 
int_t = t_func(2)- t_func(1); 


15 


del_f = f_func(2)- f_func(1) ; 
% Calculate the many "bins" to track tonal of interest. 


span = (max_v/c_air) * peak_func(1,3) ; 
span = ceil (span/del_f); 


% Locate peaks. Search for the strongest tonal in the vicinity of 
% tonal's frequency and time. ‘time intervals' before and after 
% tonal's time is as input in 'no_point'. 


a = zeros (1, (length(no_point)*2 + 1)); 
% matrix to store freq index of tonal. 


a(no_point + 1) = peak_func (1,5); % tonal's freq_index. 
% locate tonals before and after tO 
fori = 1: no_point 

[b,n] = max (lofar_func( (peak_func(1,4)- 


i 
peak_func( 
a(no_point + 1 - i) = peak_func(1,5) + n-1; 


aes 
1,5) : peak_func(1,5)+ span)); 


[c,m] = max (lofar_func( (peak_func(1,4)+i) , ... 
peak_func(1,5)- span : peak_func(1,5))); 
a(no_point + 1 + i) = peak_func(1,5) - span + m-1; 
end 


df_dt = diff(a) * del_f/int_t: 


% Check that df_dt is negative for all points. 
b =(df_dt>=0); % Flag if positive 


b = sum(b); 
if b == 

flag = 1; % All negative 
else 

flag = 0; % Not all negative 
end 


df_dt = mean (df_dt); 


DWV %%% 0% VV %%%%% end of function %%VWhhhV%%%%%% %o% 
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APPENDIX F. MATLAB FUNCTION CODE : peak_compare.m 


function [v_out,R_out,plot_check] = 
peak_compare(lofar,f_func,t_func,peak_func,time_duration,threshold_vec) 
% Finds velocity through use of correlation on input Power Spectrum. 





% Purpose : This function finds the velocity of best match to an input 
% power spectral through the use of correlation against templates of 


° 


% synthetic generated Doppler Shift patterns. These templates are of 
% generated based on inputs df/dt and velocity. The Correlation Rxy 
% is tabulated for an array of velocities. Graph of Rxy vs velocity 

% is plotted and the velocity corresponding to the highest Rxy is 

% output from the function. 

% Input Variables 

%  \ofar : Power Spectral Density Matrix 

%  f_func, t_func : Freq and time array corresponding to lofar 
%  peak_func : matrix containing information on time, frequency 
% and df/dt of tonal of interest. Format as follows : 
% 

%  PeakNo. Time Frequency Time_index Frequency_Index df/dt 
Ofc, = 4 rithasd ssh sie She Vateat se eee ea Le st 3 

% Output Variables 

% v_out : Velocity of best match 

% R_out : Rxy value of best match for v_out 

% Local Variables 

% — time_duration : Total length of time of signal to analyse (s) 
% n,m,span : counting indexes. 

% total : Index of correlation values 

% Rxy : Correlation values 

% Functions called : generate_template 


% Written by How Whye Keong, NPS, 11 Feb 2009%%%%%%%%%%%%% Vo 


° 


% Definition of constants 


° 


tic 

max_v = 100; % Max velocity of contacts expected (m/s) 
c_air = 343; % Speed of sound 

%%% Extract the segment of 'interest' around tonal of interest %%% 


% Translate time duration into no. of time blocks before CPA 
77 


n = floor((time_duration /(t_func(2) - t_func(1)))); 


% To ensure that the duration of the template signal to be generated 
% is exactly of the same length as that of the original signal. 


time_duration = (n*2*(t_func(2) - t_func(1))); 


% Calculate number of frequency blocks. 
% An arbitary factor of 1.3 is used to extend span to ensure that the 
% entire section of the Doppler shift pattern is captured 


span = ((max_v/c_air) * peak_func(1,3)) / (f_func(2)-f_func(1)); 
span = ceil(1.3*span); 


% Extract section of interest from main lofar matrix 


% \ndex of freq of interest 


f0_| = peak_func (1,5) 
1 % Index of time of interest 


5); 
tO_| = peak_func (1,4); 
t_sect = t_func(t0_I - n : t0_I + n); 
f_sect = f_func(f0_| - (span) : f0_I + (span)); 
lofar_sect = lofar(t0_I-n:t0_I+n, f0_I- span: f0_I + span); 








[m n] = size (lofar_sect); 

% Normalise lofar_sect. 
norm_fac = sum(lofar_sect,2); 
N = norm_fac * ones(1,n); 
lofar_sect = lofar_sect ./ N ; 


%%% Correlate lofar_sect to possible templates %%% 


vel_min = 30; 
vel_max = 100; 


% Compute velocity resolution based on freq resolution 


vel_resolution = (1/peak_func(1,3))* c_air; 
N = floor((vel_max-vel_min)/vel_ resolution); 


vel = linspace (vel_min , vel_max, N); 
% Generate the template based on df/dt and carry out correlation. 


df_dt = peak_func (1,6); 
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if (df_dt >= 0) || (peak_func (1,7)==0) ; 


disp('No Doppler Shift. df/dt inadmissible’); 
v_out = 99.999; 

R_out = 0; 

plot_check = 0; 


else 
total = zeros (1, length(vel)); 
for j = 1 : length(vel) 
v = vel(j); 


[lofar_template] = generate_template(v, time_duration , ... 
25000 ,peak_func(3),df_dt); 
lofar_template = lofar_template(: , f0_I1- span: f0_l... 


+ span); 


Rxy = xcorr2(lofar_sect, lofar_template); 
total(j)= max(max(Rxy)); 


end 


index = ((max(total)) == total); 
Rxy = total(index); 


% determine frequency band of f_0 and set threshold. 
diff_vec = threshold_vec (1,:) - peak_func(1,3); 
[m n] = min(abs(diff_vec)); 


if Rxy >= threshold_vec(4,n) 


v_out = vel(index); 
R_out = total(index); 


disp(’ '); 

disp(['Contact with Doppler Shift found att =',... 
num2str(peak_func(1,2)),'s.']); 

disp(['Velocity of contact is '\num2str(v_out),' m/s.']); 


% Produce plot 
plot_check = 2; 
lofar_db_1 = 10* log1O(lofar); 
lofar_db = 10*log10 (lofar_sect); 
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figure; 


subplot(2,2,2) 

pcolor (f_func(1:400 
t_func(peak_func 

peak_func(1,4) + 
lofar_db_1(peak_func(1,4)- 10:... 
peak_func(1,4) + 10, 1:400)); 

shading interp 

xlabel('Frequency (Hz)') 

ylabel('Time (s)’) 

colorbar 

title((LOFARgram’); 


—~woa 
Se 


subplot(2,2,3) 

pcolor (f_sect,t_sect,lofar_db) 
shading interp 
xlabel('Frequency (Hz)') 
ylabel('Time (s)') 

colorbar 

grid on; 

title("Tonal Analysed’); 


subplot(2,2,4) 
plot(vel,total); 
xlabel('Velocity (m/s)’'); 
ylabel('Rxy'); 
grid on; 
title(‘Plot of Rxy vs Velocity’); 
annotation(‘textbox','String',{'Threshold :’,... 
num2str(threshold_vec(4,n),4)}.... 
‘FitHeightToText’,'off.... 
‘Position’ ,[0.8422 0.3773 0.05937 0.05674)); 


else 
v_out = vel(index); 
R_out = total(index); 


plot_check = 1; 


disp(' ') 
disp('No Doppler Shift Detected. Rxy less than Threshold’); 


end 
end 


toc 


DWV %%%%%%%% End of Function %%%%%%%%%%%%%% 
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APPENDIX G. MATLAB FUNCTION CODE : generate_template.m 


function [lofar_out] = generate_template(vel,t_duration,Fs,f_0,df_dt) 
% This function generates template with Doppler Shift pattern based of 
% input velocity. 


% Purpose : This function generates a template of Doppler shift 


% pattern. 

% The pattern is centred on the input freq f_O and df/dt. 

% The function uses these two parameters to generate the 
% time-domain signal. The power spectrum is computed using 
% the 'time_to_freq' function. 

% Input Variables 

% vel : Velocity 

%  t_duration : length of sample to generate in time 

% f.0 : Frequency a CPA 

%  df/dt : Rate of change at CPA 

% Output Variables 

%  lofar_out: Power Spectrum with Doppler Shift pattern 

% f,t : frequency and time matrix correponding to 

% lofar_out 


% Local Variables 


% t : Time (s) 

% fF Frequency (Hz) 
% FL : Range (m) 

% ?p Pressure (pa) 


% Functions called : time_to freq 
% Written by How Whye Keong, NPS, 11 Feb 2009 %%%%%%HV%%HWVMNV%%%%%%% Vo 


c = 343; % Sound of Speed in Air 
A = 10; % Amplitude (Arbitary) 


t = linspace ((-t_duration/2),(t_duration/2),(t_duration * Fs)) ; 
% Calculate the Pressure as function of time. 
% r=vel* sqrt( (t.42)+((vel*f_0)/(c * df_dt))*2); 
% as required to check CPA range 
f = f_0 * (1 - ((vel.*t)./ (sqrt( (c*t).42 + (vel*f_O/df_dt)42)))); 


p =A* sin(2*pi.*f.*t); 
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% Convert the pressure-time profile into frequency domain 
[f,t,lofar_out] = time_to_freq(p',Fs,4,1); 


DWV VV %%%%%%%%%% End of Function %%%%%%%%%% Vo %o Vo %o%o %o%o 
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APPENDIX H. OUTPUT RESULT REPORTS FROM auto_Doppler_detector.m 


1. Run 1 — Threshold yp + 150, Ry, Threshold py, + 0.50; 


Results from analysis of 07_aug_08.txt. 
Threshold for energy based detection = 0.1093. 
No of peaks found = 131. 

No of Doppler Shifts detected = 10. 


Start time : 8/7/2008 10:20:52.55. 


Peak No. Time Freq Result Velocity Rxy Range 
1 74.71 122.07 Contact with Doppler Shift found. 48.26 56.43 106.80 
2 623.42 94.60 Contact with Doppler Shift found. 53.33 51.63 168.48 
3 814.13 125.12 No Doppler Shift. df/dt inadmissible. 
4 816.10 91.55 No Doppler Shift Detected. Rxy less than Threshold. 
5 818.07 122.07 No Doppler Shift. df/dt inadmissible. 
6 822.98 103.76 No Doppler Shift Detected. Rxy less than Threshold. 
7 835.76 115.97 No Doppler Shift Detected. Rxy less than Threshold. 
8 879.01 103.76 No Doppler Shift Detected. Rxy less than Threshold. 
9 1089.56 103.76 No Doppler Shift Detected. Rxy less than Threshold. 
10 1092.51 109.86 No Doppler Shift. df/dt inadmissible. 
11 1215.56 112.92 No Doppler Shift. df/dt inadmissible. 
12 1230.30 103.76 No Doppler Shift. df/dt inadmissible. 
13 1232.27 103.76 No Doppler Shift. df/dt inadmissible. 
14 1276.51 115.97 No Doppler Shift. df/dt inadmissible. 
15 1422.98 112.92 No Doppler Shift Detected. Rxy less than Threshold. 
16 1457.39 91.55 No Doppler Shift Detected. Rxy less than Threshold. 
17 1490.81 109.86 No Doppler Shift. df/dt inadmissible. 
18 1516.54 125.12 Contact with Doppler Shift found. 41.67 19.14 17.00 
19 1519.49 103.76 No Doppler Shift Detected. Rxy less than Threshold. 
20 1526.37 115.97 No Doppler Shift Detected. Rxy less than Threshold. 
21 1575.52 122.07 No Doppler Shift. df/dt inadmissible. 
22 1582.37 137.33 Contact with Doppler Shift found. 61.11 47.38 321.10 
23 1610.91 115.97 No Doppler Shift. df/dt inadmissible. 
24 1619.76 100.71 No Doppler Shift Detected. Rxy less than Threshold. 
25 1636.47 91.55 No Doppler Shift. df/dt inadmissible. 
26 1659.08 91.55 No Doppler Shift Detected. Rxy less than Threshold. 
27 1795.73 250.24 Contact with Doppler Shift found. 32.80 8.36 12.97 
28 1805.73 97.66 No Doppler Shift Detected. Rxy less than Threshold. 
29 1811.63 91.55 No Doppler Shift Detected. Rxy less than Threshold. 
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1859.80 
1872.58 
1874.54 
1924.68 
1931.56 
1943.35 
1946.30 
1950.24 
1957.12 
2004.30 
2028.88 
2214.85 
2119.46 
2223.69 
2229.59 
2232.54 
2263.02 
2269.90 
2305.29 
2311.18 
2371.15 
2390.81 
2467.66 
2601.35 
2611.18 
2616.10 
2632.81 
2635.76 
2642.64 
2656.40 
2692.78 
2716.54 
2726.37 
2735.22 
2739.15 
2754.88 
2790.27 
2796.17 
2822.71 
2827.63 
2849.25 
2865.96 
2873.83 
2882.68 
2920.03 
2938.71 


125.12 
125.12 
97.66 
137.33 
109.86 
97.66 
134.28 
97.66 
94.60 
125.12 
109.86 
91.55 
131.23 
125.12 
91.55 
125.12 
91.55 
94.60 
119.02 
100.71 
91.55 
97.66 
100.71 
97.66 
91.55 
91.55 
109.86 
97.66 
106.81 
97.66 
161.74 
106.81 
109.86 
125.12 
125.12 
161.74 
112.92 
91.55 
91.55 
125.12 
112.92 
125.12 
91.55 
91.55 
97.66 
106.81 


No Doppler Shift. df/dt inadmissible. 
No Doppler Shift. df/dt inadmissible. 
No Doppler Shift. df/dt inadmissible. 
Contact with Doppler Shift found. 48.15 
No Doppler Shift. df/dt inadmissible. 


No Doppler Shift Detected. Rxy less than Threshold. 


Contact with Doppler Shift found. 38.08 


No Doppler Shift Detected. Rxy less than Threshold. 


No Doppler Shift. df/dt inadmissible. 
No Doppler Shift. df/dt inadmissible. 


No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 


Contact with Doppler Shift found. 44.00 
No Doppler Shift. df/dt inadmissible. 


No Doppler Shift Detected. Rxy less than Threshold. 


No Doppler Shift. df/dt inadmissible. 
No Doppler Shift. df/dt inadmissible. 


No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 


No Doppler Shift. df/dt inadmissible. 


No Doppler Shift Detected. Rxy less than Threshold. 


No Doppler Shift. df/dt inadmissible. 


No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 


No Doppler Shift. df/dt inadmissible. 
No Doppler Shift. df/dt inadmissible. 
No Doppler Shift. df/dt inadmissible. 


No Doppler Shift Detected. Rxy less than Threshold. 


No Doppler Shift. df/dt inadmissible. 
No Doppler Shift. df/dt inadmissible. 
No Doppler Shift. df/dt inadmissible. 


No Doppler Shift Detected. Rxy less than Threshold. 


No Doppler Shift. df/dt inadmissible. 


No Doppler Shift Detected. Rxy less than Threshold. 


No Doppler Shift. df/dt inadmissible. 


No Doppler Shift Detected. Rxy less than Threshold. 


No Doppler Shift. df/dt inadmissible. 
No Doppler Shift. df/dt inadmissible. 


No Doppler Shift Detected. Rxy less than Threshold. 


No Doppler Shift. df/dt inadmissible. 


No Doppler Shift Detected. Rxy less than Threshold. 
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30.36 


18.78 


46.24 


85.42 


16.62 


119.29 


2943.62 
2948.54 
2956.40 
2971.15 
2983.93 
2987.86 
3008.68 
3011.63 
3013.59 
3023.42 
3037.19 
3048.00 
3050.95 
3056.85 
3060.78 
3069.63 
3073.56 
3075.52 
3078.47 
3083.39 
3095.19 
3105.02 
3113.86 
3121.73 
3125.66 
3204.30 
3209.22 
3259.35 
3264.27 
3278.03 
3280.98 
3286.88 
3292.78 
3296.71 
3300.81 
3303.76 
3314.58 
3317.53 
3348.98 
3381.42 
3385.35 
3389.29 
3424.68 
3443.35 
3447.29 
3451.22 


97.66 
91.55 
125.12 
91.55 
125.12 
91.55 
97.66 
125.12 
100.71 
91.55 
97.66 
112.92 
128.17 
122.07 
94.60 
91.55 
103.76 
91.55 
115.97 
100.71 
122.07 
91.55 
125.12 
97.66 
125.12 
112.92 
97.66 
125.12 
91.55 
106.81 
106.81 
97.66 
115.97 
125.12 
125.12 
91.55 
97.66 
106.81 
100.71 
97.66 
97.66 
100.71 
283.81 
97.66 
125.12 
161.74 


No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 


No Doppler Shift. df/dt inadmissible. 
No Doppler Shift. df/dt inadmissible. 
No Doppler Shift. df/dt inadmissible. 


No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 


No Doppler Shift. df/dt inadmissible. 


No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 


No Doppler Shift. df/dt inadmissible. 


No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 


No Doppler Shift. df/dt inadmissible. 


No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 


Contact with Doppler Shift found. 50.42 


No Doppler Shift Detected. Rxy less than Threshold. 


No Doppler Shift. df/dt inadmissible. 


No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 


No Doppler Shift. df/dt inadmissible. 
No Doppler Shift. df/dt inadmissible. 


No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 


No Doppler Shift. df/dt inadmissible. 
No Doppler Shift. df/dt inadmissible. 


No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 


No Doppler Shift. df/dt inadmissible. 


No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 


No Doppler Shift. df/dt inadmissible. 
No Doppler Shift. df/dt inadmissible. 
Contact with Doppler Shift found. 30.00 


No Doppler Shift Detected. Rxy less than Threshold. 


No Doppler Shift. df/dt inadmissible. 
No Doppler Shift. df/dt inadmissible. 
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18.28 


9.79 


35.14 


19.19 


122 
123 
124 
125 
126 
127 
128 
129 
130 
131 


3456.13 
3461.05 
3463.02 
3485.62 
3489.56 
3511.18 
3516.10 
3522.98 
3526.91 
3622.44 


100.71 
97.66 
125.12 
125.12 
112.92 
109.86 
109.86 
125.12 
115.97 
125.12 


No Doppler Shift. df/dt inadmissible. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift. df/dt inadmissible. 
No Doppler Shift. df/dt inadmissible. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift. df/dt inadmissible. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift. df/dt inadmissible. 


86 


2. Run 2 — Threshold yp + 150, Rxy Threshold up, + 1.00, 


Results from analysis of 07_aug_08.txt. 
Threshold for energy based detection = 0.1093. 
No of peaks found = 131. 

No of Doppler Shifts detected = 5. 


Start time : 8/7/2008 10:20:52.55. 


Peak No. Time Freq Result Velocity Rxy Range 
1 74.71. 122.07 Contact with Doppler Shift found. 48.26 56.43 106.80 
2 623.42 94.60 Contact with Doppler Shift found. 53.33 51.63 168.48 
3 814.13 125.12 No Doppler Shift. df/dt inadmissible. 
4 816.10 91.55 No Doppler Shift Detected. Rxy less than Threshold. 
5 818.07 122.07 No Doppler Shift. df/dt inadmissible. 
6 822.98 103.76 No Doppler Shift Detected. Rxy less than Threshold. 
7 835.76 115.97 No Doppler Shift Detected. Rxy less than Threshold. 
8 879.01 103.76 No Doppler Shift Detected. Rxy less than Threshold. 
9 1089.56 103.76 No Doppler Shift Detected. Rxy less than Threshold. 
10 1092.51 109.86 No Doppler Shift. df/dt inadmissible. 
11 1215.56 112.92 No Doppler Shift. df/dt inadmissible. 
12 1230.30 103.76 No Doppler Shift. df/dt inadmissible. 
13 1232.27 103.76 No Doppler Shift. df/dt inadmissible. 
14 1276.51 115.97 No Doppler Shift. df/dt inadmissible. 
15 1422.98 112.92 No Doppler Shift Detected. Rxy less than Threshold. 
16 1457.39 91.55 No Doppler Shift Detected. Rxy less than Threshold. 
17 1490.81 109.86 No Doppler Shift. df/dt inadmissible. 
18 1516.54 125.12 No Doppler Shift Detected. Rxy less than Threshold. 
19 1519.49 103.76 No Doppler Shift Detected. Rxy less than Threshold. 
20 1526.37 115.97 No Doppler Shift Detected. Rxy less than Threshold. 
21 1575.52 122.07 No Doppler Shift. df/dt inadmissible. 
22 1582.37 137.33 Contact with Doppler Shift found. 61.11 47.38 321.10 
23 1610.91 115.97 No Doppler Shift. df/dt inadmissible. 
24 1619.76 100.71 No Doppler Shift Detected. Rxy less than Threshold. 
25 1636.47 91.55 No Doppler Shift. df/dt inadmissible. 
26 1659.08 91.55 No Doppler Shift Detected. Rxy less than Threshold. 
27 1795.73 250.24 No Doppler Shift Detected. Rxy less than Threshold. 
28 1805.73 97.66 No Doppler Shift Detected. Rxy less than Threshold. 
29 1811.63 91.55 No Doppler Shift Detected. Rxy less than Threshold. 
30 1859.80 125.12 No Doppler Shift. df/dt inadmissible. 
31 1872.58 125.12 No Doppler Shift. df/dt inadmissible. 
32 1874.54 97.66 No Doppler Shift. df/dt inadmissible. 
33 1924.68 137.33 Contact with Doppler Shift found. 48.15 30.36 85.42 
34 1931.56 109.86 No Doppler Shift. df/dt inadmissible. 
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1943.35 
1946.30 
1950.24 
1957.12 
2004.30 
2028.88 
2214.85 
2119.46 
2223.69 
2229.59 
2232.54 
2263.02 
2269.90 
2305.29 
2311.18 
2371.15 
2390.81 
2467.66 
2601.35 
2611.18 
2616.10 
2632.81 
2635.76 
2642.64 
2656.40 
2692.78 
2716.54 
2726.37 
2735.22 
2739.15 
2754.88 
2790.27 
2796.17 
2822.71 
2827.63 
2849.25 
2865.96 
2873.83 
2882.68 
2920.03 
2938.71 
2943.62 
2948.54 
2956.40 
2971.15 
2983.93 


97.66 
134.28 
97.66 
94.60 
125.12 
109.86 
91.55 
131.23 
125.12 
91.55 
125.12 
91.55 
94.60 
119.02 
100.71 
91.55 
97.66 
100.71 
97.66 
91.55 
91.55 
109.86 
97.66 
106.81 
97.66 
161.74 
106.81 
109.86 
125.12 
125.12 
161.74 
112.92 
91.55 
91.55 
125.12 
112.92 
125.12 
91.55 
91.55 
97.66 
106.81 
97.66 
91.55 
125.12 
91.55 
125.12 


No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift. df/dt inadmissible. 

No Doppler Shift. df/dt inadmissible. 

No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 
Contact with Doppler Shift found. 

No Doppler Shift. df/dt inadmissible. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift. df/dt inadmissible. 

No Doppler Shift. df/dt inadmissible. 

No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift. df/dt inadmissible. 

No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift. df/dt inadmissible. 

No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift. df/dt inadmissible. 

No Doppler Shift. df/dt inadmissible. 

No Doppler Shift. df/dt inadmissible. 

No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift. df/dt inadmissible. 

No Doppler Shift. df/dt inadmissible. 

No Doppler Shift. df/dt inadmissible. 

No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift. df/dt inadmissible. 

No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift. df/dt inadmissible. 

No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift. df/dt inadmissible. 

No Doppler Shift. df/dt inadmissible. 

No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift. df/dt inadmissible. 

No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift. df/dt inadmissible. 

No Doppler Shift. df/dt inadmissible. 

No Doppler Shift. df/dt inadmissible. 
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44.00 46.24 119.29 


2987.86 
3008.68 
3011.63 
3013.59 
3023.42 
3037.19 
3048.00 
3050.95 
3056.85 
3060.78 
3069.63 
3073.56 
3075.52 
3078.47 
3083.39 
3095.19 
3105.02 
3113.86 
3116.81 
3121.73 
3125.66 
3204.30 
3209.22 
3264.27 
3278.03 
3280.98 
3286.88 
3292.78 
3296.71 
3300.81 
3303.76 
3314.58 
3317.53 
3348.98 
3381.42 
3385.35 
3389.29 
3443.35 
3447.29 
3451.22 
3456.13 
3461.05 
3463.02 
3485.62 


91.55 
97.66 
125.12 
100.71 
91.55 
97.66 
112.92 
128.17 
122.07 
94.60 
91.55 
103.76 
91.55 
115.97 
100.71 
122.07 
91.55 
125.12 
125.12 
97.66 
125.12 
112.92 
97.66 
91.55 
106.81 
106.81 
97.66 
115.97 
125.12 
125.12 
91.55 
97.66 
106.81 
100.71 
97.66 
97.66 
100.71 
97.66 
125.12 
161.74 
100.71 
97.66 
125.12 
125.12 


No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift. df/dt inadmissible. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift. df/dt inadmissible. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift. df/dt inadmissible. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift. df/dt inadmissible. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift. df/dt inadmissible. 
No Doppler Shift. df/dt inadmissible. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift. df/dt inadmissible. 
No Doppler Shift. df/dt inadmissible. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift. df/dt inadmissible. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift. df/dt inadmissible. 
No Doppler Shift. df/dt inadmissible. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift. df/dt inadmissible. 
No Doppler Shift. df/dt inadmissible. 
No Doppler Shift. df/dt inadmissible. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift. df/dt inadmissible. 
No Doppler Shift. df/dt inadmissible. 
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126 
127 
128 
129 
130 
131 


3489.56 
3511.18 
3516.10 
3522.98 
3526.91 
3622.44 


112.92 
109.86 
109.86 
125.12 
115.97 
125.12 


No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift. df/dt inadmissible. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift. df/dt inadmissible. 
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3: Run 3 — Threshold yp + 200, Rxy Threshold up, + 0.50, 


Results from analysis of 07_aug_08.txt. 


Threshold for energy based detection = 0.1380. 


No of peaks found = 113. 


No of Doppler Shifts detected = 7. 


Start time : 8/7/2008 10:20:52.55. 


Peak No. Time 
1 74.71 
2 623.42 
3 814.13 
4 816.10 
5 818.07 
6 822.98 
7 835.76 
8 1089.56 
9 1092.51 
10 1215.56 
11 1273.56 
12 1276.51 
13 1422.98 
14 1430.84 
15 1457.39 
16 1519.49 
17 1521.46 
18 1524.41 
19 1575.52 
20 1582.41 
21 1619.76 
22 1653.18 
23 1794.74 
24 1805.73 
25 1811.63 
26 1872.58 
27 1874.54 
28 1924.68 
29 1931.56 
30 1943.35 
31 1946.30 
32 1957.12 
33 2004.30 
34 2028.88 
35 2119.49 


Freq 
122.07 
94.60 
125.12 
91.55 
122.07 
103.76 
115.97 
103.76 
109.86 
112.92 
100.71 
115.97 
112.92 
106.81 
91.55 
103.76 
106.81 
91.55 
122.07 
137.33 
100.71 
94.60 
225.83 
97.66 
91.55 
125.12 
97.66 
137.33 
109.86 
97.66 
134.28 
94.60 
125.12 
109.86 
131.23 


Result Velocity Rxy Range 
Contact with Doppler Shift found. 48.26 56.43 106.80 
Contact with Doppler Shift found. 53.33 51.63 168.48 


No Doppler Shift. df/dt inadmissible. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift. df/dt inadmissible. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift. df/dt inadmissible. 
No Doppler Shift. df/dt inadmissible. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift. df/dt inadmissible. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift. df/dt inadmissible. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift. df/dt inadmissible. 


Contact with Doppler Shift found. 61.11 47.38 321.10 


No Doppler Shift Detected. Rxy less than Threshold. 

No Doppler Shift Detected. Rxy less than Threshold. 

No Doppler Shift Detected. Rxy less than Threshold. 

No Doppler Shift Detected. Rxy less than Threshold. 

No Doppler Shift Detected. Rxy less than Threshold. 

No Doppler Shift. df/dt inadmissible. 

No Doppler Shift. df/dt inadmissible. 

Contact with Doppler Shift found. 48.15 30.36 
No Doppler Shift. df/dt inadmissible. 

No Doppler Shift Detected. Rxy less than Threshold. 

Contact with Doppler Shift found. 38.08 18.78 
No Doppler Shift. df/dt inadmissible. 

No Doppler Shift. df/dt inadmissible. 

No Doppler Shift Detected. Rxy less than Threshold. 

Contact with Doppler Shift found. 44.00 46.24 
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85.42 


16.62 


119.29 


2214.85 
2216.81 
2229.59 
2232.54 
2263.02 
2269.90 
2301.35 
2305.29 
2311.18 
2371.15 
2467.66 
2601.35 
2611.18 
2616.10 
2632.81 
2635.76 
2642.64 
2656.40 
2692.78 
2716.54 
2733.25 
2735.22 
2739.15 
2790.27 
2796.17 
2812.88 
2822.71 
2825.66 
2849.25 
2865.96 
2873.83 
2882.68 
2920.03 
2938.71 
2943.62 
2948.54 
2983.93 
2987.86 
3008.68 
3013.59 
3023.42 
3041.12 
3048.00 
3050.95 
3056.85 
3060.78 


91.55 
97.66 
91.55 
125.12 
91.55 
94.60 
161.74 
119.02 
100.71 
91.55 
100.71 
97.66 
91.55 
91.55 
109.86 
97.66 
106.81 
97.66 
161.74 
106.81 
97.66 
125.12 
125.12 
112.92 
91.55 
94.60 
91.55 
115.97 
112.92 
125.12 
91.55 
91.55 
97.66 
106.81 
97.66 
91.55 
125.12 
91.55 
97.66 
100.71 
91.55 
91.55 
112.92 
128.17 
122.07 
94.60 


No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift. df/dt inadmissible. 

No Doppler Shift. df/dt inadmissible. 

No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift. df/dt inadmissible. 

No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift. df/dt inadmissible. 

No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift. df/dt inadmissible. 

No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 


No Doppler Shift. 
No Doppler Shift. 
No Doppler Shift. 
No Doppler Shift. 
No Doppler Shift. 
No Doppler Shift. 


df/dt inadmissible. 
df/dt inadmissible. 
df/dt inadmissible. 
df/dt inadmissible. 
df/dt inadmissible. 
df/dt inadmissible. 


No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift. df/dt inadmissible. 

No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift. df/dt inadmissible. 

No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift. df/dt inadmissible. 

No Doppler Shift. df/dt inadmissible. 

No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift. df/dt inadmissible. 

No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift. df/dt inadmissible. 

No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 
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3069.63 
3073.56 
3075.52 
3078.47 
3083.39 
3095.19 
3105.02 
3113.86 
3121.73 
3204.30 
3209.22 
3262.30 
3264.27 
3278.03 
3280.98 
3286.88 
3292.78 
3300.81 
3303.76 
3389.29 
3443.35 
3447.29 
3451.22 
3456.13 
3461.05 
3463.02 
3485.62 
3489.56 
3511.18 
3516.10 
3526.91 
3622.44 


91.55 
103.76 
91.55 
115.97 
100.71 
122.07 
91.55 
125.12 
97.66 
112.92 
97.66 
97.66 
91.55 
106.81 
106.81 
97.66 
115.97 
125.12 
91.55 
100.71 
97.66 
125.12 
161.74 
100.71 
97.66 
125.12 
125.12 
112.92 
109.86 
109.86 
115.97 
125.12 


No Doppler Shift. df/dt inadmissible. 


No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 


No Doppler Shift Detected. Rxy less than Threshold 
Contact with Doppler Shift found. 


50.42 18.28 35.14 


No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 


No Doppler Shift. df/dt inadmissible. 


No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 


No Doppler Shift. df/dt inadmissible. 


No Doppler Shift Detected. Rxy less than Threshold. 


No Doppler Shift. df/dt inadmissible. 


No Doppler Shift Detected. Rxy less than Threshold. 


No Doppler Shift. df/dt inadmissible. 
No Doppler Shift. df/dt inadmissible. 
No Doppler Shift. df/dt inadmissible. 


No Doppler Shift Detected. Rxy less than Threshold. 


No Doppler Shift. df/dt inadmissible. 
No Doppler Shift. df/dt inadmissible. 


No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 


No Doppler Shift. df/dt inadmissible. 
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4. Run 4 - Threshold yp + 200, Rxy Threshold up, + 1.00; 


Results from analysis of 07_aug_08.txt. 
Threshold for energy based detection = 0.1380. 
No of peaks found = 113. 

No of Doppler Shifts detected = 5. 


Start time : 8/7/2008 10:20:52.55. 


Peak No. Time Freq Result Velocity Rxy Range 
1 74.71. 122.07 Contact with Doppler Shift found. 48.26 56.43 106.80 
2 623.42 94.60 Contact with Doppler Shift found. 53.33 51.63 168.48 
3 814.13 125.12 No Doppler Shift. df/dt inadmissible. 
4 816.10 91.55 No Doppler Shift Detected. Rxy less than Threshold. 
5 818.07 122.07 No Doppler Shift. df/dt inadmissible. 
6 822.98 103.76 No Doppler Shift Detected. Rxy less than Threshold. 
7 835.76 115.97 No Doppler Shift Detected. Rxy less than Threshold. 
8 1089.56 103.76 No Doppler Shift Detected. Rxy less than Threshold. 
9 1092.51 109.86 No Doppler Shift. df/dt inadmissible. 
10 1215.56 112.92 No Doppler Shift. df/dt inadmissible. 
11 1273.56 100.71 No Doppler Shift Detected. Rxy less than Threshold. 
12 1276.51 115.97 No Doppler Shift. df/dt inadmissible. 
13 1422.98 112.92 No Doppler Shift Detected. Rxy less than Threshold. 
14 1430.84 106.81 No Doppler Shift. df/dt inadmissible. 
15 1457.39 91.55 No Doppler Shift Detected. Rxy less than Threshold. 
16 1519.49 103.76 No Doppler Shift Detected. Rxy less than Threshold. 
17 1521.46 106.81 No Doppler Shift Detected. Rxy less than Threshold. 
18 1524.41 91.55 No Doppler Shift Detected. Rxy less than Threshold. 
19 1575.52 122.07 No Doppler Shift. df/dt inadmissible. 
20 1582.41 137.33 Contact with Doppler Shift found. 61.11 47.38 321.10 
21 1619.76 100.71 No Doppler Shift Detected. Rxy less than Threshold. 
22 1653.18 94.60 No Doppler Shift Detected. Rxy less than Threshold. 
23 1794.74 225.83 No Doppler Shift Detected. Rxy less than Threshold. 
24 1805.73 97.66 No Doppler Shift Detected. Rxy less than Threshold. 
25 1811.63 91.55 No Doppler Shift Detected. Rxy less than Threshold. 
26 1872.58 125.12 No Doppler Shift. df/dt inadmissible. 
27 1874.54 97.66 No Doppler Shift. df/dt inadmissible. 
28 1924.68 137.33 Contact with Doppler Shift found. 48.15 30.36 85.42 
29 1931.56 109.86 No Doppler Shift. df/dt inadmissible. 
30 1943.35 97.66 No Doppler Shift Detected. Rxy less than Threshold. 
31 1946.30 134.28 No Doppler Shift Detected. Rxy less than Threshold. 
32 1957.12 94.60 No Doppler Shift. df/dt inadmissible. 
33 2004.30 125.12 No Doppler Shift. df/dt inadmissible. 
34 2028.88 109.86 No Doppler Shift Detected. Rxy less than Threshold. 
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2119.49 
2214.85 
2216.81 
2229.59 
2232.54 
2263.02 
2269.90 
2301.35 
2305.29 
2311.18 
2371.15 
2467.66 
2601.35 
2611.18 
2616.10 
2632.81 
2635.76 
2642.64 
2656.40 
2692.78 
2716.54 
2733.25 
2735.22 
2739.15 
2790.27 
2796.17 
2812.88 
2822.71 
2825.66 
2849.25 
2865.96 
2873.83 
2882.68 
2920.03 
2938.71 
2943.62 
2948.54 
2983.93 
2987.86 
3008.68 
3013.59 
3023.42 
3041.12 
3048.00 
3050.95 
3056.85 


131.23 
91.55 
97.66 
91.55 

125.12 
91.55 
94.60 

161.74 

119.02 

100.71 
91.55 

100.71 
97.66 
91.55 
91.55 

109.86 
97.66 

106.81 
97.66 

161.74 

106.81 
97.66 

125.12 

125.12 

112.92 
91.55 
94.60 
91.55 

115.97 

112.92 

125.12 
91.55 
91.55 
97.66 

106.81 
97.66 
91.55 

125.12 
91.55 
97.66 

100.71 
91.55 
91.55 

112.92 

128.17 

122.07 


Contact with Doppler Shift found. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift. df/dt inadmissible. 
No Doppler Shift. df/dt inadmissible. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift. df/dt inadmissible. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift. df/dt inadmissible. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift. df/dt inadmissible. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift. df/dt inadmissible. 
No Doppler Shift. df/dt inadmissible. 


44.00 46.24 119.29 


No Doppler Shift. 
No Doppler Shift. 
No Doppler Shift. 


df/dt inadmissible. 
df/dt inadmissible. 
df/dt inadmissible. 


No Doppler Shift. df/dt inadmissible. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift. df/dt inadmissible. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift. df/dt inadmissible. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift. df/dt inadmissible. 
No Doppler Shift. df/dt inadmissible. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift. df/dt inadmissible. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift. df/dt inadmissible. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 
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3060.78 
3069.63 
3073.56 
3075.52 
3078.47 
3083.39 
3095.19 
3105.02 
3113.86 
3121.73 
3204.30 
3209.22 
3262.30 
3264.27 
3278.03 
3280.98 
3286.88 
3292.78 
3300.81 
3303.76 
3389.29 
3443.35 
3447.29 
3451.22 
3456.13 
3461.05 
3463.02 
3485.62 
3489.56 
3511.18 
3516.10 
3526.91 
3622.44 


94.60 
91.55 
103.76 
91.55 
115.97 
100.71 
122.07 
91.55 
125.12 
97.66 
112.92 
97.66 
97.66 
91.55 
106.81 
106.81 
97.66 
115.97 
125.12 
91.55 
100.71 
97.66 
125.12 
161.74 
100.71 
97.66 
125.12 
125.12 
112.92 
109.86 
109.86 
115.97 
125.12 


No Doppler Shift Detected. 


Rxy less than Threshold. 


No Doppler Shift. df/dt inadmissible. 


No Doppler Shift Detected. 
No Doppler Shift Detected. 
No Doppler Shift Detected. 
No Doppler Shift Detected. 
No Doppler Shift Detected. 
No Doppler Shift Detected. 
No Doppler Shift Detected. 
No Doppler Shift Detected. 
No Doppler Shift Detected. 
No Doppler Shift Detected. 
No Doppler Shift Detected. 


Rxy less than Threshold. 
Rxy less than Threshold. 
Rxy less than Threshold. 
Rxy less than Threshold. 
Rxy less than Threshold. 
Rxy less than Threshold. 
Rxy less than Threshold. 
Rxy less than Threshold. 
Rxy less than Threshold. 
Rxy less than Threshold. 
Rxy less than Threshold. 


No Doppler Shift. df/dt inadmissible. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift. df/dt inadmissible. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift. df/dt inadmissible. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift. df/dt inadmissible. 
No Doppler Shift. df/dt inadmissible. 
No Doppler Shift. df/dt inadmissible. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift. df/dt inadmissible. 
No Doppler Shift. df/dt inadmissible. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift Detected. Rxy less than Threshold. 
No Doppler Shift. df/dt inadmissible. 
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5. Run 5 — Threshold yp + 150, Rxy Threshold u, + 0.50, 


Results from analysis of 07_aug_08.txt. 
Threshold for energy based detection = 0.3000. 
No of peaks found = 65. 

No of Doppler Shifts detected = 5. 


Start time : 8/7/2008 10:20:52.55. 


Peak No. Time Freq Result Velocity Rxy Range 
1 74.71 122.07 Contact with Doppler Shift found. 48.26 56.43 106.80 
2 623.42 94.60 Contact with Doppler Shift found. 53.33 51.63 168.48 
3 1276.51 115.97 No Doppler Shift. df/dt inadmissible. 
4 1422.98 112.92 No Doppler Shift Detected. Rxy less than Threshold. 
5 1427.90 94.60 No Doppler Shift Detected. Rxy less than Threshold. 
6 1524.41 91.55 No Doppler Shift Detected. Rxy less than Threshold. 
7 1580.44 170.90 Contact with Doppler Shift found. 46.97 41.53 177.04 
8 1582.41 149.54 Contact with Doppler Shift found. 61.11 47.38 321.10 
9 1653.18 94.60 No Doppler Shift Detected. Rxy less than Threshold. 
10 1924.68 137.33 Contact with Doppler Shift found. 48.15 30.36 85.42 
11 1927.63 112.92 No Doppler Shift Detected. Rxy less than Threshold. 
12 1929.59 91.55 No Doppler Shift Detected. Rxy less than Threshold. 
13 1931.56 109.86 No Doppler Shift. df/dt inadmissible. 
14 1933.52 109.86 No Doppler Shift. df/dt inadmissible. 
15 1935.49 100.71 No Doppler Shift Detected. Rxy less than Threshold. 
16 2214.85 91.55 No Doppler Shift Detected. Rxy less than Threshold. 
17 2216.81 97.66 No Doppler Shift Detected. Rxy less than Threshold. 
18 2229.59 91.55 No Doppler Shift Detected. Rxy less than Threshold. 
19 2259.08 100.71 No Doppler Shift Detected. Rxy less than Threshold. 
20 2263.02 91.55 No Doppler Shift. df/dt inadmissible. 
21 2269.90 94.60 No Doppler Shift Detected. Rxy less than Threshold. 
22 2301.35 161.74 No Doppler Shift. df/dt inadmissible. 
23 2305.29 119.02 No Doppler Shift Detected. Rxy less than Threshold. 
24 2311.18 100.71 No Doppler Shift Detected. Rxy less than Threshold. 
25 2467.66 100.71 No Doppler Shift Detected. Rxy less than Threshold. 
26 2601.35 97.66 No Doppler Shift. df/dt inadmissible. 
27 2611.18 91.55 No Doppler Shift Detected. Rxy less than Threshold. 
28 2616.10 91.55 No Doppler Shift. df/dt inadmissible. 
29 2635.76 97.66 No Doppler Shift Detected. Rxy less than Threshold. 
30 2642.64 106.81 No Doppler Shift Detected. Rxy less than Threshold. 
31 2716.54 106.81 No Doppler Shift. df/dt inadmissible. 
32 2790.27 112.92 No Doppler Shift Detected. Rxy less than Threshold. 
33 2793.22 106.81 No Doppler Shift Detected. Rxy less than Threshold. 
34 2796.17 91.55 No Doppler Shift. df/dt inadmissible. 


97 


2815.83 
2822.71 
2830.57 
2833.52 
2849.25 
2873.83 
2882.68 
2884.64 
2943.62 
2948.54 
2950.51 
2952.47 
2987.86 
3002.78 
3008.68 
3023.42 
3041.12 
3043.08 
3048.00 
3050.95 
3060.78 
3078.47 
3095.19 
3111.90 
3113.86 
3121.73 
3278.03 
3292.78 
3456.13 
3516.10 
3526.91 


94.60 
91.55 
106.81 
91.55 
112.92 
91.55 
91.55 
97.66 
97.66 
91.55 
106.81 
91.55 
91.55 
112.92 
97.66 
91.55 
91.55 
97.66 
112.92 
128.17 
94.60 
115.97 
122.07 
112.92 
125.12 
97.66 
106.81 
115.97 
100.71 
109.86 
115.97 


No Doppler Shift Detected. 
No Doppler Shift Detected. 
No Doppler Shift Detected. 
No Doppler Shift Detected. 
No Doppler Shift Detected. 


Rxy less than Threshold. 
Rxy less than Threshold. 
Rxy less than Threshold. 
Rxy less than Threshold. 
Rxy less than Threshold. 


No Doppler Shift. df/dt inadmissible. 


No Doppler Shift Detected. 
No Doppler Shift Detected. 
No Doppler Shift Detected. 
No Doppler Shift Detected. 


Rxy less than Threshold. 
Rxy less than Threshold. 
Rxy less than Threshold. 
Rxy less than Threshold. 


No Doppler Shift. df/dt inadmissible. 


No Doppler Shift Detected. 
No Doppler Shift Detected. 
No Doppler Shift Detected. 
No Doppler Shift Detected. 
No Doppler Shift Detected. 
No Doppler Shift Detected. 


Rxy less than Threshold. 
Rxy less than Threshold. 
Rxy less than Threshold. 
Rxy less than Threshold. 
Rxy less than Threshold. 
Rxy less than Threshold. 


No Doppler Shift. df/dt inadmissible. 


No Doppler Shift Detected. 
No Doppler Shift Detected. 
No Doppler Shift Detected. 
No Doppler Shift Detected. 
No Doppler Shift Detected. 
No Doppler Shift Detected. 
No Doppler Shift Detected. 
No Doppler Shift Detected. 
No Doppler Shift Detected. 
No Doppler Shift Detected. 


Rxy less than Threshold. 
Rxy less than Threshold. 
Rxy less than Threshold. 
Rxy less than Threshold. 
Rxy less than Threshold. 
Rxy less than Threshold. 
Rxy less than Threshold. 
Rxy less than Threshold. 
Rxy less than Threshold. 
Rxy less than Threshold. 


No Doppler Shift. df/dt inadmissible. 


No Doppler Shift Detected. 
No Doppler Shift Detected. 


Rxy less than Threshold. 
Rxy less than Threshold. 
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